AFRL-AFOSR-UK-TR-201 6-0021 


Distributed  Low  Temperature  Combustion  133024 


Peter  Lindstedt 

IMPERIAL  COLLEGE  OF  SCIENCE  TECHNOLOGY  &  MEDICINE 


09/07/2016 
Final  Report 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Air  Force  Research  Laboratory 
AF  Office  Of  Scientific  Research  (AFOSR)/  IOE 
Arlington,  Virginia  22203 
Air  Force  Materiel  Command 


FORM  SF  298 


Page  1  of  1 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Executive  Services,  Directorate  (0704-0188).  Respondents  should  be  aware  that 
notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

07-09-2016  Final  01  Feb  2013  to  31  Jul  2016 

4.  TITLE  AND  SUBTITLE 

Distributed  Low  Temperature  Combustion:  Fundamental  Understanding  of  Combustion 
Regime  Transitions 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

FA8655- 13- 1-3024 

5c.  PROGRAM  ELEMENT  NUMBER 

61102F 

6.  AUTHOR(S) 

Peter  Lindstedt 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

IMPERIAL  COLLEGE  OF  SCIENCE  TECHNOLOGY  &  MEDICINE  REPORT  NUMBER 

EXHIBITION  RD 
LONDON,  SW7  2BT  GB 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

EOARD 
Unit  4515 

APO  AE  09421-4515 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

A  DISTRIBUTION  UNLIMITED:  PB  Public  Release 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  objective  of  the  current  study  is  to  bring  fundamental  understanding  of  the  impact  of  the  chemical  (Tau_c)  and  flow  (Tau_f)  timescales  on 
combustion  regime  transitions  in  turbulent  premixed  flows.  Hence,  Tau_c  is  here  varied  via  the  mixture  stoichiometry  (Phi)  with  variations  in 
Tau_f  pursued  in  a  parallel  study.  Aerodynamically  stabilised  dimethyl  ether  (DME)  flames  in  a  backstep  burnt  opposed  jet  configuration, 
featuring  fractal  generated  multi-scale  turbulence(Re_t  >  350),  are  used  to  study  decoupled  parameters  affecting  Tau_c  and  combustion 
regime  transitions  from  conventional  flamelets  into  the  distributed  reaction  zone  regime.  The  choice  of  DME  is  partly  due  to  the  potential 
practical  relevance,  but  also  due  to  the  fundamentally  different  chemical  behaviour  as  compared  to  ethanol.  The  latter  fuel  has  also  been 
considered  along  with  methane.  Work  has  also  been  performed  on  the  further  assessment  of  chemical  mechanisms  for  the  considered  fuels 
(e.g.  DME)  to  establish  their  ability  to  reproduce  laminar  flame  and  auto-ignition  properties.  The  chemical  mechanisms  where  then  used  to 
determine  parameters  required  in  the  estimation  of  the  different  Damkoehler  numbers.  This  chemistry  can  be  used  to  identify  the  impact  of  the 
major  chemical  pathways  on  combustion  mode  transitions.  The  conceptual  multifluid  approach  of  Spalding  can  be  used  to  avoid  the 
limitations  associated  with  the  common  bimodal  (two-fluid)  description  and  is  here  explored  via  simultaneous  OH-PLIF  and  PIV,  permitting  the 
identification  of  five  separate  fluid  states. _ 

15.  SUBJECT  TERMS 

EOARD,  Low  Temperature  Combustion,  Combustion  Regime  Transitions 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

CUMMINGS,  RUSSELL 

Unclassified 

Unclassified 

Unclassified 

SAR 

75 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

01 1-44-1895-616021 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

AFRL/AFOSR  lOE 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

AFRL-AFOSR-UK-TR-201 6-0021 


Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.1 8 


https://livelink.ebs.afrl.af.mil/liveliiik/llisapi.dll 


10/13/2016 


Final  Report  FA8655-13-1-3024 


Distributed  Low  Temperature  Combustion: 
Fundamental  Understanding  of  Combustion 
Regime  Transitions 


Period  of  Performance:  1  February  2013  to  31  July  2016 

F.  Hampp,  P.  Kraus,  P.  Simatos  and  R.  P.  Lindstedt 
Department  of  Mechanical  Engineering,  Imperial  College, 
Exhibition  Road,  London  SI/1/7  2AZ,  UK 


’Principal  Investigator,  email:  p.lindstedt@imperial.ac.uk 


1 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


2 


Contents 

1  Introduction .  8 

2  Methods,  Assumptions,  and  Procedures .  9 

2.1  Multi-Fluid  Description  .  11 

2.2  Classification  of  the  Multiple  Fluid  States .  12 

2.2.1  Density  Segregation  Technique  .  12 

2.2.2  Fluid  Detection  using  OH-PLIF .  13 

2.2.3  Hot  Combustion  Products  .  13 

2.2.4  Mixing/Preheat  (Low  Temperature  Chemistry)  Fluids  .  14 

2.2.5  Strongly  Reacting  Fluid  .  14 

2.2.6  Product  Fluid .  15 

2.2.7  Mildly  Reacting  Fluid .  16 

2.3  Sensitivity  Analysis .  17 

2.4  Spatial  Multi-Fluid  Resolution .  17 

3  Results  and  Discussion .  18 

3.1  Multi-Fluid  Analysis .  18 

3.2  Unconditional  Velocity  Statistics .  20 

3.3  Conditional  Velocity  Statistics .  20 

3.3.1  Conditional  Reactant  Fluid  Velocity .  21 

3.3.2  Conditional  Mixing  Fluid  Velocity .  21 

3.3.3  Conditional  Mildly  Reacting  Fluid  Velocity .  22 

3.3.4  Conditional  Strongly  Reacting  Fluid  Velocity .  22 

3.3.5  Conditional  Strain  Distribution  on  Material  Surfaces .  22 

3.3.6  Strain  Distribution  along  the  Reactant  Fluid  Surface  .  23 

3.3.7  Strain  Distribution  along  Chemically  Active  Fluid  Surfaces  ....  23 

3.4  Bimodal  Flow  Analysis .  24 

3.5  Impact  of  Excess  Enthalpy  on  the  Burning  Mode  Transition .  25 

3.5.1  Multi-Fluid  Probability  Statistics .  25 

3.5.2  Conditional  Velocity  Statistics .  26 

3.5.3  Conditional  Strain  along  Material  Surfaces .  27 

3.6  Combustion  Regime  Transitions .  28 

3.7  Calculation  Methods  for  Distributed  Combustion  Regimes .  30 

3.7.1  Efficiency  Gains  from  Parallelisation  .  32 

4  Conclusions .  33 

5  Acknowledgements .  35 

6  References .  35 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


3 


7  Tables .  42 

8  Figures  .  46 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


4 


List  of  Figures 

1  Experimental  configuration .  47 

2  Coordinate  system  convention  and  reference  windows .  47 

3  Identification  of  the  multiple  fluid  states  in  a  sample  Mie  scattering  and 

OH-PLIF  image .  48 

4  Schematic  visualisation  of  the  multi-fluid  post-processing  methodology.  .  49 

5  OH  concentration  following  mixing  of  cold  air  with  HCP  to  establish  the 

maximum  OH  signal  levels  in  the  absence  of  fuel .  50 

6  Theoretical  mixture  seeding  density  due  to  UN  and  LN  stream  mixing  as  a 

function  of  blending  quantity  (i.e.  mixing  temperature)  and  seeding  density 
ratio  for  and  HCP  temperature  =  1700  K .  51 

7  DME  auto-ignition  delay  times .  51 

8  Laminar  flame  calculations  to  determine  the  extinction  strain  and  to  define 

the  strongly  reacting  fluid  threshold .  52 

9  Laminar  flame  calculations  to  visualise  the  change  in  burning  mode  with 

increasing  strain  in  the  back-to-burnt  configuration  for  premixed  DME/air 
flames  at  varying  $ .  53 

10  Laminar  flame  calculations  in  a  BTB  geometry  for  varying  HCP  conditions.  54 

1 1  OH  lifetime  in  the  post  combustion  environment  for  DME  at  <f>  =  0.60, 

0.80  and  1.00  under  extinction  conditions .  55 

12  Sensitivity  analysis  on  the  fluid  separating  threshold  values  for  varying  <f>.  55 

13  Example  of  a  quinary  multi-fluid  field .  56 

14  Multi-fluid  probability  statistics  for  DME  cases  at  $  =  0.2  -1.0  aligned  at 

the  first  thermal  alteration  iso-contour .  57 

15  Multi-fluid  interface  diagrams  for  DME/air  mixtures .  58 

16  Unconditional  mean  axial  velocity  and  axial  and  radial  fluctuations  along 

the  burner  centre  axis .  59 

17  Conditional  mean  axial  reactant  fluid  velocity  and  fluctuations  along  the 

stagnation  point  streamline .  60 

18  Conditional  mean  axial  mixing  fluid  velocity  and  fluctuations  along  the 

stagnation  point  streamline .  61 

19  Conditional  mean  axial  mildly  reacting  fluid  velocity  and  fluctuations  along 

the  stagnation  point  streamline .  62 

20  Conditional  mean  axial  strongly  reacting  fluid  velocity  and  fluctuation  along 

the  stagnation  point  streamline .  63 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


5 


21  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  reactant  fluid 

surface .  64 

22  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  mildly  reacting 

fluid  surface .  64 

23  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  strongly  reacting 

fluid  surface .  65 

24  Coordinate  system  convention .  65 

25  Multi-fluid  probability  for  DME  cases  at  $  =  0.6  with  HCP  temperature  = 

1600-  1800  K .  66 

26  Multi-fluid  interface  diagram  for  DME  /  air  at  $  =  0.6  and  varying  sup¬ 
porting  HCP  temperature .  67 

27  Conditional  mean  axial  mixing  fluid  velocity  and  fluctuations  for  varying 

HCP  enthalpy  along  the  stagnation  point  streamline .  68 

28  Conditional  mean  axial  mildly  reacting  fluid  velocity  and  fluctuations  for 

varying  HCP  enthalpy  along  the  stagnation  point  streamline .  69 

29  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  reactant  fluid 

iso-contour. .  70 

30  Combustion  regime  transitions  of  DME/air  mixtures  at  varying  $  visu¬ 
alised  in  a  Borghi  diagram .  70 

31  Auto-ignition  manifold  for  DME/air  mixtures .  71 

32  The  theoretical  maximum  speed-up  according  to  Amdahl’s  Law  [64].  .  .  72 

33  The  difference  in  memory  alignment  of  a  Structure  of  Arrays  and  Arrays 

of  Structures .  72 

34  Validating  the  MPI  results  using  the  serial  code .  73 

35  The  achieved  speed-up  of  PIPER  using  MPI .  73 

List  of  Tables 

1  Summary  of  utilised  OH  intensities  for  varying  fuel  and  equivalence  ratio 

to  define  the  fluid  separating  thresholds .  43 

2  Summary  of  turbulent  and  chemical  conditions  to  derive  the  turbulent  Reynolds, 
Damkohler  and  Karlovitz  numbers  for  DME/air  mixtures  at  varying  $  and 

low  rate  of  strain .  43 

3  Summary  of  utilised  OH  intensities  for  DME/air  at  $  =  0.6  and  varying 

HCP  temperature  to  define  the  fluid  separating  thresholds .  44 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


6 


4  Summary  of  turbulent  and  chemical  conditions  to  derive  the  turbulent  Reynolds, 

Damkohler  and  Karlovitz  numbers  for  DME/air  at  $  =  0.6  and  varying 
HCP  enthalpy .  44 

5  Physical  and  resolved  length  scales  in  pm.  Xprv  and  A m f  is  the  PIV  and 

multi-fluid  resolution,  respectively.  The  current  minimum  laminar  flame 
thickness  at  4>  =  1.0  is  min(<5/).  A/)  is  the  mean  scalar  dissipation  layer 
thickness  [59]  and  Xp  the  Batchelor  scale  [60] .  44 

6  Multi-fluid  interface  probabilities  for  DME/air  mixtures  by  traversing  along 

the  stagnation  point  streamline .  45 

7  Multi-fluid  interface  probabilities  for  varying  HCP  enthalpy  by  traversing 

along  the  stagnation  point  streamline .  45 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


7 


Summary 

The  objective  of  the  current  study  is  to  bring  fundamental  understanding  of  the  impact  of 
the  chemical  (rc)  and  flow  (t/)  time-scales  on  combustion  regime  transitions  in  turbulent 
premixed  flows.  Hence,  rc  is  here  varied  via  the  mixture  stoichiometry  (<f>)  with  variations 
in  Tf  pursued  in  a  parallel  study.  Aerodynamically  stabilised  dimethyl  ether  (DME)  flames 
in  a  back-to-burnt  opposed  jet  configuration,  featuring  fractal  generated  multi-scale  tur¬ 
bulence  (Ret  >  350),  are  used  to  study  decoupled  parameters  affecting  rc  and  combustion 
regime  transitions  from  conventional  flamelets  into  the  distributed  reaction  zone  regime. 
The  choice  of  DME  is  partly  due  to  the  potential  practical  relevance,  but  also  due  to  the  fun¬ 
damentally  different  chemical  behaviour  as  compared  to  ethanol.  The  latter  fuel  has  also 
been  considered  along  with  methane.  Work  has  also  been  performed  on  the  further  assess¬ 
ment  of  chemical  mechanisms  for  the  considered  fuels  (e.g.  DME)  to  establish  their  ability 
to  reproduce  laminar  flame  and  auto-ignition  properties.  The  chemical  mechanisms  where 
then  used  to  determine  parameters  required  in  the  estimation  of  the  different  Damkohler 
numbers.  This  chemistry  can  be  used  to  identify  the  impact  of  the  major  chemical  pathways 
on  combustion  mode  transitions.  The  conceptual  multi-fluid  approach  of  Spalding  can  be 
used  to  avoid  the  limitations  associated  with  the  common  bimodal  (two-fluid)  description 
and  is  here  explored  via  simultaneous  OH-PLIF  and  PIV,  permitting  the  identification  of 
five  separate  fluid  states.  The  latter  are  defined  as  reactant,  combustion  product,  mixing, 
mildly  reacting  and  strongly  reacting  (e.g.  flamelet)  fluids.  To  distinguish  between  the  dif¬ 
ferent  chemically  active  fluids,  thresholds  are  defined  based  on  the  detected  OH  signal. 
The  sensitivity  of  the  distribution  between  the  different  fluid  states  to  the  threshold  values 
is  also  evaluated.  In  order  to  correlate  the  fluid  state  thresholds  to  their  physical  inter¬ 
pretation,  laminar  flame  calculations  in  the  directly  corresponding  geometry  were  carried 
out.  The  work  suggests  that  a  consistent  treatment  of  applied  threshold  levels  is  necessary 
when  comparing  different  types  of  simulations  (e.g.  DNS)  and  experiments  for  reacting 
fluids  with  multiple  states.  However,  the  work  also  suggests  that  such  characterisations 
are  possible  and  useful  in  delineating  combustion  regime  transitions.  The  application  of 
calculation  methods  that  do  not  rely  on  a  presumed  (e.g.  flamelet)  structure  requires  the 
introduction  of  multiple  chemical  time-scales.  To  this  effect,  the  scope  of  parallelising  the 
treatment  of  the  chemistry  in  the  context  of  transported  joint-pd/f  methods  has  also  been  ex¬ 
plored  on  different  platforms  in  order  to  permit  effective  computations  featuring  extended 
mechanisms  for  practical  fuels.  The  latter  are  of  a  size  approaching  that  of  mechanisms 
based  on  the  rapid  decomposition  hypothesis  leading  to  the  C1-C2  and/or  C1-C4  founda¬ 
tion  chemistry.  It  is  shown  that  close  to  the  theoretical  maximum  speed-up  can  be  achieved. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


1  Introduction 


1  Introduction 

Combustion  technologies  operating  in  the  absence  of  distinct  flame  fronts  are  expected  to 
lead  to  improved  control,  reduced  fuel  consumption  through  stable  fuel  lean  operation,  low 
maximum  temperatures  and  reduced  emissions.  This  development  induces  the  need  to  ex¬ 
tend  laboratory  research  to  significantly  higher  turbulence  levels  and  reduced  Damkohler 
numbers  to  provide  fundamental  understanding  of  combustion  regime  transitions  including 
distributed  and/or  flameless  regimes,  e.g.  flameless  oxidation  (FLOX)  [1,  2],  high  temper¬ 
ature  air  combustion  (HiTAC)  [3]  or  moderate  and  intense  low  oxygen  dilution  (MILD) 
combustion  [4,  5],  The  opposed  jet  geometry  has  been  extensively  utilised  to  investigate 
turbulent  (and  laminar)  non-pre mixed,  partially  premixed  and  premixed  combustion  as 
well  as,  more  recently,  regime  transitions  [6,  7,  8,  9].  The  geometry  entails  decisive  advan¬ 
tages  rendering  it  indispensable  for  experimental  combustion  studies  as  well  as  computa¬ 
tional  model  development  and  validation  [10,  11].  These  are:  (i)  Excellent  optical  access 
for  laser-based  diagnostic  measurements;  (ii)  Accurate  experimental  control  of  boundary 
conditions;  (iii)  Aerodynamic  flame  stabilisation,  rather  than  through  pilot  flames,  leading 
to  flame  dynamics  and  extinction  being  directly  related  to  the  intrinsic  aerothermochem- 
istry  of  the  combustion  process;  (iv)  Individual  control  of  variables  affecting  the  chemi¬ 
cal  and  turbulent  time-scale  and  (v)  a  compact  physical  domain  ideal  for  LES  and  other 
computational  studies.  The  last  three  points  identify  the  current  experimental  facility  as 
arguably  ideal  for  investigations  of  combustion  regime  transitions. 

Dimethyl  ether  (DME)  has  been  promoted  since  the  mid-1990s  [12,  13]  as  an  attractive 
alternative  to  Diesel.  Numerous  investigations  are  found  on  the  utilisation  of  DME  in  en¬ 
gines,  e.g.  [14,  15],  yet  fundamental  flame  studies  are  relatively  rare.  Existing  studies  have 
investigated  properties  such  as  laminar  burning  velocities  [16,  17,  18]  and  extinction  strain 
rates  [18].  Recently,  Fuest  et  al.  [19]  presented  an  experimental  study  of  laminar  and  turbu¬ 
lent  DME  jet  flames  using  Raman/Rayleigh  scattering  and  CO-LIF  supported  by  laminar 
flame  calculations.  Species  and  temperature  profiles  were  obtained  for  two  equivalence 
ratios  and  good  agreement  was  found  between  experiment  and  computation.  Furthermore, 
a  higher  impact  of  differential  diffusion  was  established  for  all  DME/air  flames  compared 
to  CHi i/air  flames  due  to  significantly  increased  level  of  molecular  hydrogen.  Gabet  et 
al.  [20]  compared  the  flame  structure  of  non-premixed  DME  and  CH4  /  air  turbulent  jet 
flames.  Despite  the  increasing  number  of  publications,  no  research  effort  with  respect  to 
the  combustion  regime  transition  of  premixed  DME/air  flames  has  been  found  in  litera¬ 
ture.  The  latter  issue  is  partially  addressed  in  the  current  study  by  considering  data  that 
may  permit  an  extension  of  arguably  the  most  successful  statistical  description  of  turbulent 
premixed  combustion,  i.e.  the  two-fluid  BML  model  [21],  The  work  performed  on  DME 
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is  here  used  as  an  example  of  the  investigations  performed.  In  addition,  methane  [22,  23] 
and  ethanol  [23]  have  also  been  considered. 

The  statistical  description  of  combustion  approaching  distributed  reaction  zones  ne¬ 
cessitates  a  more  complex  methodology  as  significant  reaction  zone  broadening  and  in¬ 
termediate  fluid  states  are  expected.  Consequently,  the  current  work  utilises  a  back-to- 
burnt  (BTB)  opposed  jet  configuration  to  investigate  the  impact  of  the  chemical  timescale 
and  excess  enthalpy  on  the  combustion  regime  transition  of  lean  premixed  DME  flames. 
The  conceptual  multi-fluid  approach  of  Spalding  [25]  is  explored  via  simultaneous  OH 
-  PLIF  and  PIV  permitting  the  identification  of  various  combustion  gas  states.  A  range 
of  Damkohler  numbers  (Da)  from  the  conventional  propagating  flamelet  regime  well  into 
the  distributed  reaction  zone  regime  are  investigated.  The  latter  requires  the  application  of 
calculation  methods  that  do  not  rely  on  a  presumed  (e.g.  flamelet)  structure  and  inevitably 
leads  to  the  need  to  introduce  multiple  chemical  time-scales.  The  scope  of  parallelising  the 
treatment  of  the  chemistry  in  the  context  of  transported  joint-pc/f  methods  has  accordingly 
been  explored  in  order  to  permit  effective  computations  featuring  extended  mechanisms  for 
practical  fuels.  The  latter  are  of  a  size  approaching  that  of  mechanisms  based  on  the  rapid 
decomposition  hypothesis  leading  to  the  C1-C2  and/or  C1-C4  foundation  chemistry  [26], 

2  Methods,  Assumptions,  and  Procedures 

Combustion  regime  transitions  from  self-sustained  flames  (e.g.  flamelets)  to  distributed 
reactions  can  be  achieved  by  augmenting  the  turbulence  intensity  or  by  slowing  down  the 
combustion  chemistry.  In  the  present  work,  burning  mode  transitions  were  systematically 
enforced  by  altering  the  chemical  timescale  via  the  mixture  stoichiometry  (0.0  <  $  < 
1.0),  while  maintaining  tj  constant  with  Re*  >  350  (Re  ~  18,400).  This  corresponds 
to  conditions  beyond  the  intense  turbulence  regime  urms/ Sl  —  19  [27]  with  the  present 
values  covering  the  range  3.3  <  urms/ Sl  <  43.  Zhou  et  al.  [28]  identified  distributed 
reactions  in  stoichiometric  vitiated  CH4  jet  flames  for  Ret  >  240,  which  reduced  to  Ret  > 
130  for  $  =  0.4.  The  impact  of  the  excess  enthalpy  of  the  supporting  HCP  stream  on 
the  regime  transition  is  also  investigated.  Based  on  a  conventional  combustion  regime 
diagram  [29,  30]  the  conditions  cover  a  transition  from  (close  to)  corrugated  flamelets  to 
distributed  reaction  zones. 

The  multi-fluid  concept  of  Spalding  [25]  is  explored  via  simultaneous  Mie  scattering, 
particle  image  velocimetry  (PIV)  and  hydroxyl  planar  laser  induced  fluorescence  (OH- 
PLIF)  in  order  to  extend  bimodal  (two-fluid)  approaches  [21].  This  novel  analysis  permits 
the  statistical  evaluation  of  intermediate  fluid  states  that  are  found  in  low  Damkohler  num¬ 
ber  flows.  The  current  implementation  allows  segregations  of  five  different  fluid  states 
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(reactants,  combustion  products,  mixing  fluid,  mildly  and  strongly  reacting  fluids),  with  a 
thermochemical  interpretation  provided.  The  derived  multi-fluid  analysis  is  here  used  to 
delineate  systematic  combustion  regime  transitions  from  close  to  the  corrugated  flamelet 
regime  well  into  a  distributed  reaction  zone  regime  of  premixed  DME  by  means  of  (a) 
multi-fluid  probability  and  (b)  interface  statistics,  (c)  unconditional  and  (d)  conditional 
velocity  statistics,  (e)  underlying  strain  statistics  conditioned  upon  the  material  surfaces 
and  (f)  limitations  of  bimodal  segregation  techniques.  The  multi-fluid  analysis  is  further 
used  to  delineate  (g)  the  impact  of  the  excess  enthalpy  of  the  stabilising  hot  combustion 
product  on  the  burning  mode  transition  with  the  impact  of  the  fuel  chemistry  explored  by 
Hampp  [22]  and  strain  (bulk  and  turbulent)  by  Hampp  and  Lindstedt  [24],  The  conven¬ 
tional  combustion  regime  classification  is  extended  by  (h)  a  tentative  delineation  into  a  3D 
regime  manifold  to  accommodate  conditions  of  distributed  burning  modes. 

The  opposed  jet  burner  (Fig.  1),  originally  developed  by  Geyer  et  al.  [11],  was  ex¬ 
tensively  used  by  Geipel  et  al.  [10]  and  Goh  et  al.  [8,  31,  32,  33]  and  is  identical  to  the 
the  configuration  used  by  Hampp  et  al.  [22,  23].  A  cross  fractal  grid  (CFG;  [34])  was 
installed  50  mm  upstream  the  UN  exit,  featuring  a  BR  of  65  %  with  maximum  and  min¬ 
imum  bar  widths  of  2.0  mm  and  0.50  mm.  The  CFG  provides  increased  multiscale  tur¬ 
bulence  [8,  10,  31,  32],  Premixed  DME  /  air  mixtures  with  U b,UN  =  11.0  m/s  at  T(  -  y 
=  320  K  were  injected  through  the  upper  nozzle  (UN)  at  varying  $  =  0.20,  0.40,  0.60, 
0.80  and  1.0.  The  velocity  fluctuations  (urms  =  1.58  m/s)  and  the  integral  length  scale 
of  turbulence  (L/  =  4.1  mm)  were  measured  using  hot  wire  anemometry,  resulting  in  a 
turbulent  Reynolds  number  of  Ret  «  370  (with  kinematic  viscosity  (v)  of  1 7.9- 1  0”  6  for 
air  at  320  K). 

Premixed  Ho  ('I'  =  1.0)  flames  were  stabilised  on  a  stabilising  mesh  (FSM),  located 
100  mm  upstream  the  LN  exit,  providing  thermally  equilibrated  hot  combustion  products 
(HCP)  to  the  stagnation  plane.  A  flashback  arrestor  (FBA)  is  installed  for  safety.  The  HCP 
temperature  was  controlled  by  the  CO2  dilution  rate  within  the  range  of  1600  to  1800  K. 
The  nozzle  exit  temperature  and  its  peak-to-peak  variation  of  ~  1  %  were  measured  us¬ 
ing  an  R-type  thermocouple  with  an  estimated  time  constant  of  10  ms.  The  aerodynamic 
stabilisation  of  the  stagnation  plane  in  the  proximity  of  the  burner  centre  (mid-point  be¬ 
tween  the  nozzles)  was  realised  by  matching  the  momentum  of  the  LN  HCP  stream  and 
the  UN  premixed  fuel-air  mixture.  Thus,  the  target  Ret  necessitated  a  LN  bulk  velocity  of 
U b.LN  ~  4.3  m/s  (at  293  K).  The  coordinate  system  convention  is  shown  in  Fig.  2. 

Simultaneous  OH  -  PLIF  and  PIV  were  carried  out  using  the  BalNO.jU  crystal  based 
Raman  laser  of  Kerl  et  al.  [35]  (two  light  sheets;  281.7  nm  and  532  nm;  height  =  ID; 
thickness  <  0.5  mm;  At  =  25  its).  The  OH  -  PLIF  images  were  obtained  from  the  first 
laser  pulse.  The  interrogation  regions  for  the  OH  -  PLIF  and  PIV  cameras  were  set  to 
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51.9x38.6  mm  and  34.5x25.6  mm  (1376x  1023  pixels).  Aluminium  oxide  (AI2O3)  pow¬ 
der  ( pa  =  3900  kg/m3),  with  particle  diameters  of  da  50  =  0.44  pm  and  da  go  =  1.66  p m, 
was  used  as  velocity  tracer  particles.  Adaptive  PIV  with  decreasing  interrogation  window 
sizes  (128  x  128  to  48  x  48  with  a  75  %  overlap)  offers  the  highest  accuracy,  spatial  res¬ 
olution  and  robustness  at  the  penalty  of  significantly  increased  calculation  time  [36].  The 
determined  velocity  field  is  resolved  by  115  x  88  vectors,  providing  a  spatial  resolution 
of  0.30  mm.  The  experimental  and  diagnostic  setup  were  discussed  in  detail  by  Hampp  et 
al.  [22,  24], 

2.1  Multi-Fluid  Description 

Combustion  with  Da  <  1  can  lead  to  a  broadening  of  reaction  zones  as  observed  in  vitiated 
jet  flames  [37,  38,  28,  39]  and  DNS  [40,  41],  Burning  modes  with  distributed  reactions 
show  a  strong  degeneration  of  the  laminar  species  /  reaction  progress  correlation  [41]  and 
a  bimodal  two-fluid  description  (reactants  and  products)  with  a  negligible  probability  of 
encountering  chemically  active  states  becomes  problematic  [42].  Spalding  [25]  suggested 
a  multi-fluid  approach  that  permits  the  identification  of  various  intermediate  fluid  states. 
The  concept  was  explored  by  Hampp  [22]  and  Hampp  and  Lindstedt  [23,  24]  using  simul¬ 
taneous  Mie  scattering,  PIV  and  OH  -  PLIF  combined  with  a  purpose  written  algorithm 
(summarised  see  below)  that  detects  four  iso-contours  in  each  instantaneous  image  pair 
in  order  to  distinguish  between  five  different  fluid  states.  The  existence  of  relatively  rare 
events  required  3000  independent  realisations  to  maintain  statistical  accuracy. 

The  experimentally  observed  OH  signal  intensity  bands  were  linked  to  theoretically 
determined  OH  concentrations  to  provide  physical  interpretation.  Three  types  of  calcula¬ 
tions  were  performed  with  the  chemical  reaction  mechanism  of  Park  [43], 

1.  Freely  propagating  laminar  flames  calculations  replicated  the  experimental  condi¬ 
tions  of  the  lower  in-nozzle  flow  and  provided  a  value  for  the  OH  radical  concen¬ 
tration  at  the  LN  exit  that  is  used  as  reference  ([OH%).  The  measured  nozzle  exit 
temperature  was  matched  via  a  radiation  heat  loss  correction  [44]  of  7.2  -  8.9  %. 
The  computational  domain  was  resolved  by  660  distributed  cells  with  local  refine¬ 
ment  of  the  reaction  zone. 

2.  One-dimensional  strained  laminar  flame  calculations  [44]  were  performed  in  the 
twin  flame  mode  to  determine  the  flame  extinction  point  and  in  the  BTB  mode 
corresponding  to  the  experimental  configuration  to  provide  a  consistent  set  of  Da 
for  all  conditions.  The  computational  domain  was  resolved  by  390  distributed  cells 
with  local  refinement  of  the  CH  peak  leading  to  a  resolution  <  12  pm. 
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3.  Zero-dimensional  shock  tube  calculations  were  carried  out  to  determine  auto-ignition 
delay  times  and  perfectly  stirred  reactor  (PSR)  calculations  to  estimate  OH  lifetimes 
in  different  fluids.  The  time  step  of  the  PSR  calculation  was  adjusted  to  resolve  the 
ignition  onset  and  was  as  low  as  1  x  10”8  s  for  high  To. 

The  linking  of  experimental  data,  e.g.  OH  -  PLIF  signal  levels,  to  theoretically  determined 
radical  concentrations  via  laminar  flame  calculations,  to  support  post-processing  strategies 
becomes  essential  when  using  the  signal  intensity  of  one  reaction  zone  marker,  e.g.  OH  - 
PLIF  signal,  to  highlight  differences  within  the  intensity  of  the  chemical  activity. 

2.2  Classification  of  the  Multiple  Fluid  States 

The  selected  thermochemical  fluid  states  accessible  with  the  used  simultaneous  Mie  scat¬ 
tering,  PIV  and  OH  -  PLIF  technique  are  defined  below  and  sample  images  with  overlaid 
PIV  vectors  and  detected  iso-contours  are  shown  in  Fig.  3.  The  overall  flow  chart  to  deter¬ 
mine  the  different  fluid  states  is  depicted  in  Fig.  4. 

Reactant  Fluid:  Fresh  reactants  emerging  from  the  UN  that  have  not  undergone  thermal 
alteration  (i.e.  oxidation  or  mixing  processes)  are  detected  via  density  segregation. 

Mixing  Fluid:  A  fluid  state  without  detectable  OH  signal  that  has  been  exposed  to  a  ther¬ 
mal  change  (i.e.  mixing  of  HCP  with  reactants). 

Mildly  Reacting  Fluid:  A  fluid  state  with  modest  levels  of  OH  such  as  the  chemical  ac¬ 
tivity  of  ultra  lean  flames  that  are  sustained  by  the  thermal  support  provided  by  an 
external  enthalpy  source. 

Strongly  Reacting  Fluid:  Regions  with  a  strong  OH  signal  intensity  indicate  a  large  OH 
radical  pool  that  originates  from  self-sustained  (e.g.  flamelet)  burning  and  conven¬ 
tional  aerothermochemistry  conditions  and  extinction  criteria  apply  [45], 

Product  Fluid:  Burnt  combustion  products  (BCP)  from  the  UN  mixture  are  combined 
with  HCP  from  the  LN. 

2.2.1  Density  Segregation  Technique 

PIV  tracer  particle  based  density  segregation  (DS)  techniques  are  widely  used,  e.g.  [31, 46, 
47],  and  the  current  methodology  is  a  variant  that  enables  the  detection  of  more  general 
flow  structures.  The  DS  iso-contour  is  determined  based  on  a  binary  image  via  Moore- 
Neighbor  tracing  with  a  Jacobs  stopping  criteria  [48],  The  algorithm  permits  the  detec¬ 
tion  of  multiple,  fragmented  and  disconnected  zones.  Mixing  of  the  opposing  streams 
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and  subsequent  heat  release  can  result  in  ambiguous  determinations  and  constraints  were 
imposed  on  the  seeding  densities  (Nsj).  To  achieve  unambiguous  signal  identification 
Nsd,UN  >  1-75  -Ngd'LN  was  found  sufficient  to  ensure  a  unique  density  segregation  iso¬ 
contour.  Over-seeded  images  were  rejected  as  heat  release  associated  changes  of  Nsd,uN 
could  not  be  detected  due  to  pixel  saturation.  The  overall  rejection  rate  was  <5%. 


2.2.2  Fluid  Detection  using  OH-PLIF 

The  experimentally  determined  OH  signal  intensities  were  linked  to  theoretical  OH  radical 
concentrations  in  order  to  segregate  different  reactive  fluid  states.  Alternative  methods  are 
possible  by  using  other  chemical  species  (e.g.  CH  or  CH2O).  However,  with  increasing 
reaction  zone  broadening,  the  CH20  and  OH  layers  are  expected  to  penetrate  further  into 
the  zone  of  maximum  heat  release  (CH  peak)  [38],  potentially  improving  the  OH  marker 
capabilities.  Further,  the  current  procedure  has  the  benefit  of  simplicity  and  a  sustained  OH 
concentration  ([OH*]r)  in  the  HCP  was  used  as  reference  level  at  a  fixed  location  (dashed 
LN  rectangle  in  Fig.  2).  The  peak-to-peak  fluctuations  of  the  LN  exit  temperature,  i.e. 
Thcp  ±  1  %,  translates  to  an  uncertainty  in  the  OH  concentration  of  ~  7  %.  The  theoreti¬ 
cally  determined  [OH*]t  concentrations  for  different  Thcp  (7.38  -  10.8  x  10-3  mol/m3) 
were  combined  with  the  corresponding  experimental  OH  signal  intensity  /*  (normalised  to 
unity)  to  delineate  the  other  fluid  state  (FS)  intensities  ( Ifs,x )  v' a  Eq.  (1).  For  the  current 
purpose  a  simple  linear  correlation  between  the  fluorescence  signal  and  OH  concentration 
was  found  valid  over  a  comparatively  wide  range  (1200  <  T  (K)  <  2200)  [22]  with  the 
used  signal  intensities  listed  in  Table  1 . 


[OH  ]fs,x  Jt 

[OH* 


(1) 


2.2.3  Hot  Combustion  Products 

The  determined  reference  state  ([OH*]t)  is  well  defined.  However,  an  experimentally 
observed  OH  surplus  (I  >  /*)  was  found  in  the  mixing  layer  of  an  isothermal  (<f>  =  0.0) 
case.  The  cause  can  readily  be  analysed  by  considering  the  mixing  of  HCP  with  air  from 
the  UN  through  PSR  calculations  featuring  a  residence  time  range  covering  the  integral 
(77)  and  Kolmogorov  (tv)  timescales.  The  OH  radical  is  stable  at  the  HCP  conditions  as 
shown  by  the  normalised  OH  concentration  in  the  far  right  of  Fig.  5  for  Thcp  =  1600  - 
1800  K.  However,  a  decrease  in  temperature  by  means  of  cold  air  admixture  results  in  a 
surging  OH  concentration,  peaking  at  Tq/Thcp  ~  0.9,  followed  by  a  decay.  The  excess 
OH  concentration  ([OH]  /  [OH]*  >  1)  at  low  air  blending  ratios  with  (e.g.  Tq/Thcp  > 
0.8)  results  from  the  oxidation  of  the  residual  H2  and  CO  in  the  HCP  via  02  from  the 
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air  stream  [22],  At  higher  blending  ratios  the  OH  signal  is  increasingly  quenched  and 
vanishes  around  7 ’0 /'/ V/c /'  ~  0.55.  The  maximum  increase  in  the  OH  concentration  due 
to  the  blending  of  HCP  gas  and  air  is  around  185  %  and  constant  for  all  Thcp ■  Using 
Eq.  (1),  [OH] hcp,t  (13.7  -  20.0  x  10  3  mol/m3)  can  be  transformed  into  an  expected 
OH  signal  intensity  with  Ihcp  =  1.857*,  which  is  rounded  to  2.07"  and  subsequently 
used  in  Section  2.2.6  to  define  the  product  fluid  threshold. 

2.2.4  Mixing/Preheat  (Low  Temperature  Chemistry)  Fluids 

Regions  with  a  drop  in  seeding  density  of  the  UN  stream  (i.e.  a  thermal  alteration  of  the 
fluid),  but  without  a  detected  OH  signal  (i.e.  no  distinct  chemical  activity)  were  denoted 
as  mixing  fluid.  The  change  in  the  seeding  density  (Naj)  is  a  consequence  of  the  mixing 
of  the  UN  reactant  stream  (high  seeding  density)  with  the  LN  hot  combustion  products  as 
shown  for  Thcp  =  1700  K  in  Fig.  6.  The  seeding  density  ratio  was  varied  between  1.75 
<  Nsd,uN /Nsci,ln  <  100.  The  lower  limit  follows  from  the  minimum  seeding  density 
ratio  requirement  while  the  upper  limit  is  set  by  the  ratio  of  saturation  to  background  signal 
of  the  camera.  An  estimated  HCP  blending  fraction  of  5  ±  2  %  by  volume,  depending  on 
the  seeding  density  ratio  (see  Fig.  6),  with  cold  reactants  is  sufficient  to  detect  the  thermal 
alteration.  An  estimated  cold  gas  addition  of  40  ±  10  %  is  effective  to  quench  of  the  OH 
radical  concentration  in  the  HCP.  Refinements  in  the  upper  limit  definition  are  possible 
if  there  is  a  desire  to  identify  regions  of  low  temperature  chemistry  (e.g.  characterised  by 
CH2O).  Computations  of  auto-ignition  delay  times,  defined  by  the  steepest  temperature 
gradient  (Fig.  7),  suggest  that  a  residence  time  of  one  integral  timescale  corresponds  to  an 
auto-ignition  temperature  Tign  k,  1200  K  that  is  relatively  independent  of  <f>. 

2.2.5  Strongly  Reacting  Fluid 

At  high  Damkohler  numbers,  vigorously  burning  flames  detach  from  the  stagnation  plane 
and  conventional  aerothermochemical  conditions  apply.  To  provide  a  physical  interpreta¬ 
tion,  theoretical  OH  concentrations  for  the  strongly  reacting  fluid  were  correlated  with  a 
self-sustained  burning  mode  for  each  equivalence  ratio  via  laminar  flame  calculations  [44]. 

Chemical  Timescale  Variation  The  peak  OH  concentration  at  the  extinction  point  of 
strained  counterflow  twin  flames  provides  a  reference  point  for  the  lowest  expected  OH 
concentration  of  self-sustained  burning  ([OH]str  x)  [22].  The  OH  peak  at  extinction  (ag  $) 
is  used  to  mark  a  lower  limiting  OH  concentration  ([OH]str>x)  for  each  $  as  shown  in 
Fig.  8.  This  specifies  an  approximate  threshold  (A oH,str,x)’  above  which  self-sustained 
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burning  in  absence  of  HCP  support  is  possible.  The  determined  OH  concentration  under 
extinction  condition  is  related  to  the  reference  level  ([01T:]t)  by  Eq.  (2). 


L  J  1 

Istr,x  it  A OH,str,x  ' 


(2) 


The  observed  rates  of  strain,  discussed  in  Section  3.3.5,  prevent  self-sustained  binning  for 
mixtures  below  the  lower  flammability  limit  and  thus  the  definition  of  [OH]stI%x  is  not 
meaningful  for  Dai  <  1  (Da!  =  ( Lj  ■  Sl)  ■  ( 6f  ■  Mrms)-1  =  tj  ■  r"1).  Calculations  were 
accordingly  also  performed  in  a  back-to-burnt  (rather  than  twin-flame)  configuration  to 
infer  the  corresponding  laminar  burning  velocity  and  flame  thickness  (see  Figs.  9)  values 
required  to  estimate  the  conventional  flame  property  based  Dai.  The  extinction  strain  rates 
( aq )  can  be  compared  with  the  estimated  turbulent  strain  a/  «  3950  s-1  and  the  determined 
normal  an  and  tangential  at  strain  along  the  material  surfaces  (see  Section  3.3.5).  Hence, 
it  is  expected  that  the  HCP  will  exert  some  influence  on  all  cases.  Threshold  values  are 
summarised  in  Table  1  and  extinction  conditions  are  listed  in  Table  2. 


HCP  Temperature  Variation  The  impact  of  the  Thcp  on  the  burning  mode  is  investi¬ 
gated  for  the  case  marking  the  transition  from  the  thin  reaction  zone  to  distributed  reactions, 
i.e.  <f>  =  0.6,  with  an  adiabatic  flame  temperature  of  Tad  ~  1750  K,  Sl  =  0.21  m  s^1  and  5f 
=  0.46  mm.  The  equilibrium  OH  concentration  in  the  HCP  increases  from  7.38  x  10”3  to 
10.8  x  10  3  mol/m3  with  Thcp-  The  threshold  value  are  altered  accordingly  as  defined  in 
Eq.  (2)  and  listed  in  Table  3.  The  effect  of  Thcp  on  corresponding  laminar  flames  is  shown 
in  Fig.  10  by  means  of  BTB  calculations  at  varying  rates  of  strain.  For  high  HCP  tempera¬ 
tures,  S f  reduces  steadily  by  30  %  with  increasing  strain  compared  to  the  unstrained  value 
as  depicted  in  Fig.  10a.  By  contrast,  at  low  Thcp  (=  1600  K)  6j  initially  increases  for  a  < 
1000  s_1  by  ~20  %  and  subsequently  reduces  due  to  the  increase  in  bulk  strain.  For  high 
Thcp ,  the  SL  (Fig.  10b)  and  [OH]peafc/[OHt]T  (Fig.  10c)  reduce  until  a  <  2500  s"1  and 
are  maintained  constant  for  higher  rates  of  strain.  At  reduced  Thcp  both  values  plummet 
for  a  <  1000  s  1  and  subsequently  remain  constant.  The  Dai  is  maintained  constant  at 
approximately  1.2  for  high  Thcp-  but  reduces  significantly  to  Dai  —  0.3  for  low  Thcp  as 
shown  in  Fig.  lOd.  The  resulting  global  conditions  are  summarised  in  Table  4. 

2.2.6  Product  Fluid 

The  product  fluid  threshold  is  based  on  the  maximum  OH  concentration  in  the  HCP  and 
the  burnt  combustion  products  (BCP)  resulting  from  combustion  of  the  UN  reactants.  The 
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OH  concentration  of  the  HCP  ([OH]ifpp  j)  is  reliant  upon  the  Thcp >  see  above,  while  the 
OH  concentration  in  the  BCP  ([OH]  bcp,®)  varies  with  the  equivalence  ratio  and  residence 
time.  The  lifetime  of  the  OH  radical  in  the  BCP  was  estimated  using  perfectly  stirred 
reactor  simulations,  see  Fig.  11,  with  the  residence  time  varied  between  the  estimated 
Kolmogorov  ( tv )  and  the  integral  timescales. 

The  OH  concentration  decays  to  levels  <  [ OW\hcp,t  (i.e.  2.0-  [OH^jp)  within  3.0-tv, 
(the  multi-fluid  resolution  see  Sec.  2.4)  for  DME  at  $  =  0.60,  see  Fig.  11.  For  all  mixtures, 
the  OH  concentrations  decay  to  values  close  to  [OH]pcp  within  one  integral  timescale. 
The  comparatively  rapid  decay  suggests  that  a  simple  analysis  approach  is  adequate  for 
the  current  purposes.  Accordingly,  values  corresponding  to  the  small  scale  mixing  pro¬ 
cesses  (3.0-tv,)  were  used  to  estimate  BCP  signal  intensities  in  comparatively  fresh  prod¬ 
ucts  as  summarised  in  Tables  1  and  3.  The  combined  value  was  selected  as  [OH ]prod,x  = 
max([OH]pfCPT,  [OH]pcp$).  The  product  fluid  threshold  ( ^-OH,prod,\ )  and  (he  corre¬ 
sponding  OH  signal  intensity  ( Iprod,\ )  were  defined  via  Eq.  (3). 


A 


_  [OH]proc, 

OH’prod’x  "  [OH*]t 


Iprod,x  ^  A OH,prod,x  '  I 


(3) 


2.2.7  Mildly  Reacting  Fluid 

Regions  with  modest  OH  signal  levels  can  stem  from  (i)  ignition  events,  (ii)  the  decay  of 
OH  concentration  in  combustion  products  or  (iii)  chemically  active  material  that  is  diluted 
and  sustained  by  HCP.  Such  regions  are  not  further  delineated  and  are  denoted  as  mildly 
reacting  fluid.  The  OH  concentrations  at  the  extinction  point  ([OH]stI.jX)  for  twin  flames 
were  found  to  define  the  onset  of  a  characteristic  switch  in  the  BTB  burning  mode  and 
justify  the  implementation  of  the  mildly  reacting  fluid  into  the  multi-fluid  analysis.  The 
change  is  evident  for  the  laminar  flame  thickness  and  laminar  burning  velocity  as  shown 
in  Fig.  9  and  defines  a  characteristic  strain  (aCj$).  At  higher  rates  of  strain  (a  >  ac  <j>),  the 
maximum  OH  concentration,  Sf  and  Sp  collapse  and  become  insensitive  to  the  reactant 
mixture  composition.  This  implies  that  the  flame  is  governed  by  the  HCP  and  the  OH 
concentration  and  experimental  intensities  for  the  mildly  reacting  fluid  state  are  defined  as: 

A OH, prod, x  ^  Imild,x  ^  A OH,str,x  (A) 

The  threshold  values  and  intensity  levels  are  summarised  in  Tables  1  and  3  and  the  sensi¬ 
tivity  of  the  probability  distribution  to  the  threshold  values  is  discussed  below. 
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2.3  Sensitivity  Analysis 

The  impact  of  thresholds  on  statistics  was  explored  using  a  sensitivity  analysis  for  the  case 
close  to  a  transition  from  the  thin  reaction  zone  regime  to  distributed  reaction  (<f>  >  0.60). 
Potentially,  such  cases  are  the  most  difficult  to  quantify.  The  spatial  coordinate  for  the  data 
is  normalised  by  the  integral  length  scale  of  turbulence  and  aligned  on  the  iso-contour  xs 
that  corresponds  to  the  first  thermal  change  noted  in  the  Mie  scattering  data. 

The  product  fluid  threshold  AoH,prod,.<s>=o.6  separating  combustion  products  from  the 
mildly  reacting  fluid  was  varied  between  1.6  -  2.4.  The  range  starts  below  the  uncertainty 
(~  1.8)  associated  with  oxidation  of  the  residual  HCP  products.  The  overall  variation  is 
similar  to  the  change  (~  30  %)  in  the  OH  concentration  of  the  twin  <f>  =  0.80  flame  moving 
from  near  unstrained  conditions  (a  =  75  s-1)  to  the  extinction  point  (aq  =  2000  s-1).  Thus 
the  threshold  variation  is  substantial.  Results  shown  in  Fig.  12a  highlight  a  reduction  of 
the  peak  probability  of  the  mildly  reacting  fluid,  yet  the  spatial  extent  and  general  trend  of 
the  distribution  is  preserved.  A  value  AoH,prod,,<s><o.e  =  2.0  is  preferred  as  the  lower  bound 
(see  Section  2.2.6). 

The  sensitivity  of  the  strongly  reacting  (flamelet)  burning  mode  probability  was  inves¬ 
tigated  by  applying  a  threshold  variation  from  below  the  extinction  OH  concentration  in 
the  twin  flame  configuration  with  a  symmetric  shift  around  the  defined  thresholds,  i.e.  3.0 
<  AoH,str,<s>=o.6  <  3.8.  Translating  the  corresponding  OH  variation  into  the  BTB  con¬ 
figuration,  the  inferred  strain  rate  range  includes  the  region  of  characteristic  change  in  the 
laminar  flame  thickness  and  laminar  burning  velocity  depicted  in  Fig.  9.  The  results  of  the 
sensitivity  analysis  are  shown  in  Fig.  12b.  A  reduction  in  peak  probability  is  noted,  while 
the  spatial  extent  and  distribution  trend  remains.  The  sensitivity  can  be  expected  due  to  the 
steep  decay  of  the  OH  concentration  around  the  extinction  point.  However,  it  is  apparent 
that  even  with  a  large  variation,  the  probability  of  finding  chemically  active  fluid  supported 
by  HCP  is  around  36  ±  11  %  for  the  transitional  case  (<1>  =  0.6).  The  corresponding 
probability  of  finding  strongly  reacting  (flamelet  fluid)  is  19  ±  7  %. 

2.4  Spatial  Multi-Fluid  Resolution 

The  spatial  (planar)  resolution  of  the  multi-fluid  algorithm  was  determined  by  means  of  an 
US AF- 1951  test  target  [50].  Due  to  the  difference  in  applied  smoothing  filter  width  (i.e. 
16  pixels  for  the  Mie  scattering  image  versus  4  pixels  for  the  OH  image)  the  spatial  reso¬ 
lution  of  the  multi-fluid  analysis  is  limited  by  the  detection  algorithm  of  the  Mie  scattering 
iso-contour.  Thus,  an  image  of  the  resolution  target  was  acquired  with  the  Mie  scattering 
camera  and  same  post-processing  /  smoothing  steps  applied.  The  spatial  (planar)  resolu¬ 
tion  of  the  multi-fluid  algorithm  was  determined  to  0.25  mm,  which  corresponds  to  ~  3 
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Kolmogorov  length  scales.  Consequently,  individual  and  instantaneous  fluid  states  thinner 
than  3-Lv  are  disregarded.  An  overview  of  the  spatial  scales  is  provided  in  Table  5. 

3  Results  and  Discussion 

An  example  of  a  resulting  quinary  multi-fluid  field,  corresponding  to  the  image  pair  in 
Fig.  3,  is  depicted  in  Fig.  13.  The  multi-fluid  analysis  is  used  to  delineate  combustion 
regime  transitions  of  premixed  DME/air  flames  with  the  chemical  timescale  varied  by 
means  of  the  mixture  stoichiometry.  The  turbulence  intensity  is  maintained  constant  at 
Ref  ~  380  with  mixture  specific  details  listed  in  Table  2.  The  change  in  burning  mode  from 
distributed  reactions  towards  the  corrugated  flamelet  regime  is  first  evaluated  by  means  of 
the  (i)  fluid  state  distribution  and  (ii)  interface  statistics  with  a  Damkohler  number  (Dai) 
variation,  followed  by  turbulence-chemistry  interaction  through  (iii)  unconditional  and  (iv) 
conditional  velocity  statistics  as  well  as  (v)  conditional  velocity  gradient  statistics  along  the 
material  surfaces.  Subsequently,  the  impact  of  the  HCP  enthalpy  on  the  binning  mode  tran¬ 
sition  is  investigated.  Last,  a  (vii)  tentative  classification  of  the  burning  mode  transition  is 
proposed  in  a  3D  combustion  regime  manifold. 

3.1  Multi-Fluid  Analysis 

The  evolution  and  redistribution  of  the  selected  fluid  states  is  discussed  along  the  burner 
centre  axis,  i.e.  stagnation  point  streamline  (SPS),  as  a  function  of  Dax.  The  results  high¬ 
light  the  benefit  of  extending  bimodal  segregation  techniques  and  include  intermediate 
fluid  states,  in  particular  for  low  Da  flows.  Minor  inconsistencies  of  the  stagnation  plane 
location  may  arise  due  to  jet  momentum  matching.  Hence,  the  spatial  multi-fluid  prob¬ 
abilities  were  aligned  at  the  first  thermal  alteration  of  the  fluid,  i.e.  xs  =  0,  identified  by 
the  density  segregation  iso-contour.  Inherently,  the  reactant  fluid  probability  drops  sharply 
from  unity  to  zero  at  the  origin  as  shown  in  Fig.  14a.  However,  it  re-emerges  downstream, 
extending  to  one  integral  length  scale  of  turbulence  with  its  peak  at  ~  Lj/4.  The  effect  is 
independent  of  Dai  and  accordingly  related  to  turbulent  transport.  The  recurrence  can  be 
caused  by  large  eddies  tearing  out  pockets  of  unburnt  reactants  and/or  a  three-dimensional 
effects,  yet  remains  <  5  %  for  all  cases. 

The  probability  of  finding  mixing  fluid  is  shown  in  Fig.  14b.  A  sharp  rise  at  the  ori¬ 
gin  provides  evidence  of  the  importance  of  this  fluid  state  -  particularly  adjacent  to  the 
reactant  fluid  for  flows  with  Dai  <  TO.  At  high  reactivity,  i.e.  $  =  1.0,  the  peak  prob¬ 
ability  of  the  mixing  fluid  in  direct  proximity  of  the  origin  is  significantly  reduced.  This 
suggests  an  immediate  onset  of  chemical  activity  adjacent  to  the  reactant  fluid,  without  the 
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necessity  of  HCP  support.  Moreover,  with  gradually  increasing  mixture  reactivity  the  mix¬ 
ing  fluid  probability  is  reduced  away  from  the  origin  in  favour  of  chemical  reactions,  i.e. 
smaller  quantities  of  HCP  are  required  to  initialise  OH  producing  chemical  activity.  The 
results  further  show  that  the  spatial  extent  of  the  mixing  fluid  is  limited  by  Lj,  suggesting 
a  correlation  with  turbulent  transport. 

The  probabilities  of  encountering  regions  with  mildly  and  strongly  reacting  fluids  were 
also  estimated.  With  a  decreasing  chemical  timescale  the  mildly  reacting  fluid  (i.e.  a  fluid 
with  Imiid,<s>  <  A-OH,str,< t>)  gains  significance  as  shown  in  Fig.  14c  with  the  magnitude 
and  spatial  extent  enhanced.  The  reduction  in  peak  probability  for  mixtures  with  Da!  > 
1.0  is  a  consequence  of  the  distinct  augmentation  of  the  strongly  reacting  (flamelet)  fluid 
probability,  see  Fig.  14d.  However,  due  to  the  comparatively  high  turbulence  level,  some 
impact  of  HCP  support  remains.  The  peak  probability  of  the  strongly  reacting  fluid  reaches 
72  %  for  $  =  1.00  and  is  3.5  times  higher  than  the  corresponding  mildly  reacting  fluid.  For 
mixtures  with  $  >  0.60,  the  chemically  active  fluid  states  extend  spatially  beyond  one 
integral  length  scale  due  to  dilatation. 

The  above  findings  motivated  a  statistical  analysis  of  fluid  state  interfaces  encountered 
by  traversing  through  the  quinary  multi-fluid  field  along  the  burner  axis.  The  analysis 
provides  insights  into  the  multi-fluid  structure  and  highlights  the  need  for  thermal  support 
for  low  Da  flows.  The  diagrams  in  Fig.  15  illustrate  the  multi-fluid  interface  probabilities 
for  $  =  0.2,  0.6  and  1.0.  For  $  =  0.2  the  mixing  fluid  is  primarily  (98  %)  adjacent  to  the 
reactant  fluid.  A  distinct  change  of  the  fluid  states  adjacent  to  the  reactants  is  evident  for  $ 
=  1.0  (Dai  =  5.1)  with  the  mixing  fluid  reduced  from  98  %  to  16  %  in  favour  of  the  mildly 
reacting  (6  %),  strongly  reacting  (42  %)  and  product  (36  %)  fluids.  The  detachment  of  the 
mildly  reacting  fluid  from  the  reactants  by  mixing  and/or  product  fluid  interlayers  indicates 
the  need  for  thermal  support  to  initialise  OH  producing  chemical  activity  for  Dai  <  1.0. 
The  balance  of  fluid  states  sharing  an  interface  with  the  mildly  reacting  fluid  is  strongly 
dependent  on  the  Dai  number  as  shown  in  Table  6. 

The  observed  rates  of  strain,  discussed  below,  prevent  self-sustained  burning  for  the 
leaner  mixtures  and  thus  the  strongly  reacting  fluid  is  only  defined  for  <f>  >  0.6  (Dai  > 
1.0).  The  probability  of  the  strongly  reacting  fluid  being  adjacent  to  reactants  decreases 
from  42  %  for  $  =  1.0  (Dai  =  5.1)  to  3  %  for  $  =  0.6  (Dai  =  1-2)  as  also  shown  in  Table  6 
by  means  of  interface  probabilities.  The  interface  of  the  mixing  and  strongly  reacting 
fluid  shows  the  reversed  trend  and  dominates  in  the  transitional  case  with  Dai  —  1-0.  The 
residual  strongly  reacting  fluid  shares  interfaces  with  the  mildly  reacting  or  product  fluids. 
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3.2  Unconditional  Velocity  Statistics 

The  unconditional  velocity  statistics  along  the  stagnation  point  streamline,  i.e.  normalised 
mean  axial  ( UlUb )  velocity  and  its  fluctuations,  show  the  impact  of  the  combustion  heat  re¬ 
lease  on  the  turbulent  flow  field  (  see  Fig.  16).  The  location  of  the  UN  (LN)  is  at  X/D  =  0.50 
(-0.50).  In  the  proximity  of  the  UN  exit  (X/D  >  0.2)  the  U/Ub  ratio  is  not  affected  by  the 
mixture  reactivity  as  shown  in  the  top  row  of  Fig.  16.  The  impact  of  combustion  on  U/Ut, 
becomes  evident  at  X/D  <  0.2  by  an  eased  and  lagged  deceleration  of  the  mean  flow  with 
increasing  <l>.  This  is  caused  by  an  earlier  onset  of  chemical  reaction  and  a  more  pro¬ 
nounced  flow  acceleration  with  increased  heat  release. 

The  axial  and  radial  fluctuations  ( y/u'u ' /Ub  and  y/v’v'/Ub)  are  depicted  in  the  mid¬ 
dle  and  bottom  row  of  Fig.  16.  At  X/D  >  0.2  both  components  are  independent  of  the 
mixture  reactivity.  Further  downstream  (X/D  <  0.2)  the  velocity  fluctuations  are  reduced 
with  increasing  mixture  reactivity  with  the  peak  shifted  towards  the  LN.  A  double  peak 
of  \Zu'u' /Ub  emerges  at  <J>  >  0.80,  see  middle  row  of  Fig.  16,  which  is  not  observed  at 
lower  reactivities.  The  spatial  location  of  the  peak  closer  to  the  UN  corresponds  to  the 
mean  reaction  onset.  The  second  peak  indicates  the  location  of  the  mean  interaction  of  the 
opposing  streams.  The  strong  dilatation  effects  and  associated  flow  acceleration,  pushes 
the  stagnation  plane  towards  the  LN  as  evident  in  Fig.  16  (top  and  middle  row).  Similar 
trends  were  observed  by  Goh  et  al.  [8,  33]  at  lower  turbulence  levels.  The  data  show  that 
the  combustion  mode  transition  is  sufficiently  strong  to  directly  impact  the  flow  field. 

The  velocity  field  along  the  upper  nozzle  exit,  required  for  boundary  conditions  of 
computational  studies,  was  provided  by  Hampp  [22], 

_  1  N 

Uk,FS,i,j  =  Tr  CFS,n,i,j  *  Uk:n,i,j  V 

M  n—1 

1  N  _  2 

( u'u')k,FS,i,j  =  77  E  CFS,n,i,j  ’  (Uk,n,i,j  ~  Uk,FS,i,j )  V  i,j 

A  n=i  (5) 

l  N 

CFS,i,j  =  T7  E]  CFS,n,i,j  V  i,j 
A  n—1 


3.3  Conditional  Velocity  Statistics 

Conditional  velocity  statistics  allow  the  delineation  of  the  interactions  of  each  fluid  state 
with  the  turbulent  flow  field  as  well  providing  a  link  to  established  bimodal  segregation 
techniques,  e.g.  [31,  32,  46,  47,  51],  The  multi-fluid  classification  permits  conditioning 
for  each  fluid  state,  as  shown  in  Eq.  (5),  and  thus  allows  the  quantification  of  the  contribu- 
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tion  of  each  individual  (intermediate)  fluid  state.  In  Eq.  (5),  k  is  the  axial  or  radial  velocity 
component,  n  the  instantaneous  image,  N  the  total  number  of  images  (i.e.  3000)  and  i  and 
j  the  index  variables.  The  instantaneous  conditioning  variable  (c FS,n)  is  defined  as  unity 
within  the  individual  fluid  state  and  nil  elsewhere.  The  resulting  fluid  state  progress  vari¬ 
able  (C fs)  is  compiled  identical  to  a  conventional  bimodal  reaction  progress  variable  [21] 
and,  consequently,  the  criteria  YIfs^fs  =  1  is  maintained.  The  discussion  is  based  on 
the  aligned  profiles  and  the  data  in  physical  space  are  included  to  ease  the  validation  of  nu¬ 
merical  models.  A  minimum  of  75  vectors  was  used  for  the  determination  of  conditional 
velocities. 

3.3.1  Conditional  Reactant  Fluid  Velocity 

The  conditional  mean  axial  reactant  velocity  ( Ur/Ub )  and  its  axial  (y/( u'u')r/Ub )  and  ra¬ 
dial  (\Z{v'v')r/Ub)  fluctuations  are  shown  in  Fig.  17.  The  stronger  compression  of  the 
mean  axial  flow  with  increasing  Dai  due  to  the  advanced  reaction  onset  and  elevated  heat 
release,  observed  in  Fig.  16,  remains  (see  top  row).  The  conditional  reactant  fluid  velocity 
fluctuations  (see  middle  and  bottom  rows)  show  a  significant  increase  with  decreasing  $ 
due  to  the  shift  of  the  reaction  onset  towards  the  stagnation  plane  and  thus  growing  impact 
of  the  opposing  stream. 

3.3.2  Conditional  Mixing  Fluid  Velocity 

The  mixing  fluid  probability  distribution  (Fig.  14b)  has  highlighted  the  importance  of 
this  fluid  state,  in  particular  for  combustion  processes  with  Dai  <  1-0.  The  conditional 
mean  axial  mixing  fluid  velocity  UmiX/Ub  and  turbulent  fluctuations  (^/ (u'u')miX/Ub, 
\J vr ) mix / U b)  are  depicted  in  Fig.  18.  The  mixing  velocities  are  significantly  higher 
than  the  reactant  velocities  (compare  to  Fig.  17)  for  Dai  <  1-  which  suggests  HCP  ad¬ 
dition  leading  to  increased  momentum  in  the  direction  towards  the  UN.  For  flows  with 
Dai  >  1.0,  a  distinct  drop  in  UmiX/Ub  to  the  level  of  Ur/Ub  is  evident.  Thus,  the  mixing 
fluid  is  aligned  with  the  reactant  flow  field  for  $  =  1.0,  indicating  that  only  small  HCP 
quantities  are  blended  with  the  reactant  fluid  prior  combustion. 

At  low  Dai,  the  reactants  accommodate  more  heat  addition  (HCP  blending)  while  the 
turbulence  intensity  suppresses  the  immediate  onset  of  OH  producing  chemical  reactions. 
Hence,  increased  amounts  of  combustion  products  that  alter  or  govern  the  mixing  fluid 
flow  dynamics  lead  to  a  gradual  alignment  of  UmiX/Ub  with  the  HCP  fluid  flow  direction. 
This  highlights  the  increasing  need  for  thermal  support  to  initialise  the  oxidation  process 
for  flows  with  Dai  <1.0  and  suggests  a  gradual  combustion  regime  transition  away  from 
self-sustained  (e.g.  corrugated  flamelet)  burning  towards  a  supported  distributed  mode. 
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The  corresponding  axial  (y/ (u'u')miX/Ub)  and  radial  (y/ ( v'v')miX/Ub )  fluctuations  (see 
middle  and  bottom  rows)  show  a  reduction  in  turbulent  fluctuations  with  increasing  <l>. 

3.3.3  Conditional  Mildly  Reacting  Fluid  Velocity 

Results  obtained  using  velocity  conditioning  on  the  mildly  reacting  fluid,  i.e.  chemically 
active  regions  with  reduced  OH  signal  intensity,  are  shown  in  Fig.  19.  Similarly  to  the 
mixing  fluid,  the  conditional  mean  axial  velocity  (top  row)  reveals  an  increasingly  HCP 
disturbed/driven  flow  with  decreasing  Dai.  The  impact  is  illustrated  by  the  wide  separation 
and  increasingly  positive  Umud/Ub.  The  mixture  with  $  =  1.0  leads  to  negative  Umud/Ub 
close  to  the  UN,  i.e.  is  in  line  with  the  reactant  flow.  This  implies  that  chemically  active 
regions  arise  with  little  thermal  support  from  the  HCP  stream.  For  leaner  mixtures  increas¬ 
ing  HCP  addition  disturbs/drives  the  flow  with  Umud/Ub  becoming  quickly  positive  and 
directed  towards  the  UN. 

The  conditional  velocity  fluctuations  are  shown  in  the  middle  and  bottom  rows  of 
Fig.  19.  The  axial  component  (y/ {u'u')mnd/Ub)  shows  a  distinctly  decreasing  trend  with 
increasing  $  and  exhibits  a  more  moderate  slope  with  increasing  reactivity  that  can  be 
attributed  to  dilatation.  The  radial  fluctuations  show  no  distinct  trend.  The  observations 
highlight  the  reduced  influence  of  HCP  addition  with  increasing  mixture  reactivity. 

3.3.4  Conditional  Strongly  Reacting  Fluid  Velocity 

The  conditional  strongly  reacting  fluid  velocity  and  turbulent  fluctuations,  see  Fig.  20, 
were  evaluated  for  Dai  >1.0  (i.e.  <!>  =  0.60  and  1.0).  A  clear  flow  acceleration  towards  the 
stagnation  plane  of  Ustr/Ub  is  observed,  i.e.  more  negative  values  compared  to  Umud/Ub , 
due  to  the  enhance  dilatation.  While  the  mildly  reacting  fluid  velocity  was  governed  by 
the  HCP  stream  at  X/D  <  0,  the  strongly  reacting  fluid  velocity  of  the  stoichiometric  case 
is  aligned  with  the  reactant  flow  direction  (i.e.  negative  velocities).  The  conditional  axial 
fluctuations  of  the  strongly  reacting  fluid,  see  middle  row  of  Fig.  20,  are  distinctly  reduced 
for  Dai  >  1.0  as  a  consequence  of  the  stronger  heat  release.  Cases  with  Dai  >  1  fur¬ 
ther  exhibit  slightly  lower  velocity  fluctuations  in  the  strongly  reacting  fluid  than  in  the 
corresponding  mildly  reacting  fluid,  suggesting  an  increasingly  dilatation  driven  flow. 

3.3.5  Conditional  Strain  Distribution  on  Material  Surfaces 

The  velocity  gradients  were  conditioned  upon  the  fluid  state  material  surfaces  (/?).  The 
strain  rate  (e^  =  0.5(3 Ui/dxj  +  duj/dxi ))  and  vorticity  (w.y  =  dui/dxj  —  duj/dxi) 
tensors  were  determined  from  the  instantaneous  planar  PIV  fields  and  subsequently  condi¬ 
tioned  upon  /3:  ep^j  =  e-ij  \/3  and  u)ptij  =  u>ij\/3,  where  /3  =  1  at  the  iso-contour  trajectory 
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and  nil  elsewhere.  The  tangent  to  the  iso-contour  was  estimated  at  all  positions  (()  by  a  3rd 
order  polynomial  that  was  evaluated  in  the  range  of  ±  16  pixels  along  3  (twice  the  vector 
spacing).  The  strain  rate  tensor  was  rotated  by  the  angle  0  (=  7r/2  -  6)  via  fp  =  R  •  e^, 
where  R  is  the  rotation  matrix  and  9  is  the  angle  between  the  tangential  on  jB\C  and  the 
SPS.  The  normal  (an)  and  tangential  (at)  strain  rate  components  as  well  as  the  total  strain 
(ad)  are  defined  in  Eq.  (6). 


=  7/3,11 

at  =  f/3,22  (6) 

ad  =  e^n  +  e/3,22 

Results  presented  here  include  the  movement  of  the  stagnation  point  [10]  with  an  applied 
radial  limit  of  ±  0.5  Lj  away  from  the  stagnation  point  streamline.  A  more  detailed  dis¬ 
cussion  has  been  presented  by  Hampp  and  Lindstedt  [24]. 

3.3.6  Strain  Distribution  along  the  Reactant  Fluid  Surface 

The  probability  density  function  (PDF)  of  the  aM  along  the  reactant  fluid  surface,  depicted 
in  Fig.  21,  decreases  progressively  from  mean  compressive  values  found  in  the  stagnation 
plane  proximity  (—1550  s-1)  for  Dai  <  1-0  to  below  the  bulk  strain  rate  based  on  the 
nozzle  separation  (—440  s-1)  for  Dai  =5.1.  The  PDF(a„)  is  skewed  towards  extensive 
strains  with  its  mode  shifted  from  —1680  to  —330  s-1  with  increasing  <1>.  This  can  be 
attributed  to  the  rapid  transition  to  chemically  active  fluids  states  and  a  reaction  onset 
shifted  away  from  the  stagnation  plane.  The  augmented  detachment  with  increasing  Dai 
causes  a  modest  reduction  of  the  mean  tangential  component  (950  <  at  <  690  s-1)  with 
a  distinctly  squeezed  PDF  spread  (i.e.  rms  reduced  by  ~  30  %)  and  reversed  mode  shift 
(550  -  690  s-1).  The  zero  crossings,  elevated  skewness  and  reduced  spread  of  the  PDF(a„) 
and  PDF(cit)  with  increasing  Dai  has  been  observed  in  experiments  [52]  and  DNS  [53]. 
The  dilatation  driven  flow  acceleration  for  Dai  >  TO  is  highlighted  by  the  total  rate  of 
strain  shown  in  Fig.  21  (bottom  left)  with  a  shift  from  negative  to  positive  mean  values 
(—750  <  ad  <  310  s_1).  The  earlier  reaction  onset  further  results  in  a  reduced  mean 
vorticity  of  the  order  of  50  %  for  Dai  =  5.1  compared  to  Dai  <1-0  (bottom  right  of 
Fig.  21).  The  mixing  fluid  is  found  primarily  adjacent  to  the  reactant  fluid  with  qualitatively 
the  same  strain  rate  trends. 

3.3.7  Strain  Distribution  along  Chemically  Active  Fluid  Surfaces 

For  Dai  <  1.0,  the  mildly  reacting  fluid  was  spatially  separated  (>  90  %)  from  the  reactant 
fluid  by  a  mixing  and/or  product  fluid  interlayer  that  suggests  a  need  for  thermal  support. 
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The  strong  impact  of  the  HCP  on  the  mildly  reacting  fluid  surface  results  in  a  Dai  inde¬ 
pendent  an  ~  —1340  s_1,  see  Fig.  22.  The  determined  an  is  well  beyond  the  extinction 
strain  rate  for  the  transitional  case  with  Dai  —  1.0  ( aq  =  600  s-1).  The  tails  of  the  PDF(a„) 
also  exceed  the  extinction  point  for  mixtures  with  Dai  =  5.1  (aq  =  3100  s_1).  The  lack  of 
flow  alignment  for  low  Dai  flows  was  also  found  by  Chakraborty  and  Swaminathan  [53] 
at  Re(  ~  50.  The  turbulent  mixing  driven  flow  and  thermal  HCP  support  for  the  mildly 
reacting  fluid  results  in  vorticity  levels  found  in  the  proximity  of  the  stagnation  plane  that 
are  independent  of  Dai  ■ 

The  transition  to  strongly  exothermic  fluid  states  collocates  the  flow  acceleration  nor¬ 
mal  to  the  flame  front  and  results  in  a  compressive  strain  reduction  from  —1050  to  —660  s_1, 
i.e.  around  a  40  %  reduction  for  $  =  1.0  (Dai  =  5.1)  compared  to  <!>  =  0.6  (Dai  =  1.2), 
as  depicted  in  Fig.  23.  This  further  corresponds  to  a  decrease  of  50  %  compared  to  the 
MR  surface  for  <!>  =  1.0.  The  tangential  acceleration  caused  by  combustion  results  in 
increased  extensive  strain  (650  <  at  <  830  s-1).  The  shift  of  the  PDf •’(«„)  towards  re¬ 
duced  normal  strain  with  a  similar  mean,  but  reduced  spread  of  the  PDF(aj)  coupled  with 
a  more  pronounced  skewness  with  increasing  Dai,  was  qualitatively  also  shown  in  DNS 
data  by  Chakraborty  and  Swaminathan  [53],  A  transition  from  contracting  to  dilating  ad 
(—400  to  160  s-1)  is  evident  for  high  Dai  as  shown  in  the  bottom  left  of  Fig.  23.  The 
preferential  alignment  of  the  flame  front  normal  perpendicular  to  the  extensive  strain  [54] 
becomes  evident  at  high  Dai  due  to  the  favoured  stabilisation  of  self-sustained  flames  in 
low  compressive  strain  regions  (an  <  aq  )  that  are  detached  from  the  stagnation  plane.  The 
reduction  in  turbulent  flame  speed  for  lower  Dai  [32]  leads  to  a  flame  anchoring  in  the 
proximity  of  the  stagnation  plane  with  elevated  normal  strain  levels  and  attenuated  dilata¬ 
tion  [55].  The  more  ordered  flow  field  with  increasing  Dai  is  also  reflected  by  a  significant 
vorticity  mitigation  (bottom  right  in  Fig.  23)  of  the  order  of  25  %  for  $  =  1.0  (Dai  =  5.1) 
compared  to  $  =  0.6  (Dai  =  1-2).  This  also  corresponds  to  a  ~40  %  reduction  of  the  mean 
strongly  reacting  fluid  vorticity  compared  to  the  mildly  reacting  value  for  $  =  1.0. 

3.4  Bimodal  Flow  Analysis 

Bimodal  velocity  statistics  are  frequently  inferred  from  PIV  particle  density  segregation 
techniques  (DST)  [46]  or  based  on  the  OH-PLIF  signal  [56].  Hampp  and  Lindstedt  [24] 
applied  a  PIV  based  segregation  technique  to  delineate  the  current  data  series.  A  distinctly 
reduced  gradient  transport  was  shown  with  increasing  <i>,  yet  the  transition  to  counter¬ 
gradient  transport  was  suppressed  even  for  the  stoichiometric  flame.  The  accuracy  of  two- 
fluid  descriptions  can  be  evaluated  by  comparing  the  PIV  and  OH-PLIF  based  techniques 
and  computing  the  differences  (A)  of  the  conditional  reactant  (AUr),  product  (A Uv)  and 
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slip  (A Us)  velocities  and  the  scalar  flux  (A (cu,))  as  shown  in  Fig.  24.  At  high  Da,  the 
definitions  show  good  agreement  as  only  thin  interfaces  separate  reactants  from  OH  rich 
combustion  products.  At  Da  <  1  substantial  discrepancies  emerge  as  the  Mie  scattering 
iso-contour  detaches  from  the  OH  producing  chemically  active  zones  due  to  the  need  for 
thermal  support. 

3.5  Impact  of  Excess  Enthalpy  on  the  Burning  Mode 
Transition 

The  analysis  above  has  shown  that  the  HOP  exerts  a  distinct  impact  on  the  burning  mode 
and  thus  the  impact  of  the  excess  enthalpy  is  investigated.  Local  flame  quenching  in  the 
BTB  configuration  is  primarily  caused  by  the  high  strain  acting  on  the  flame  front.  Re¬ 
ignition  occurs  either  due  to  conventional  flame  propagation  or  by  means  of  HCP  support 
(e.g.  resulting  in  auto-ignition).  The  impact  of  the  HCP  temperature  is  investigated  for  the 
case  with  Dai  —  1.0  (DME  at  $  =  0.6),  which  marks  the  approximate  transition  from  thin 
flamelets  to  distributed  reaction  zones.  The  adiabatic  flame  temperature  of  the  laminar,  un¬ 
strained,  flame  is  Tad  =  1755  K,  but  reduces  quickly  to  1555  K  under  extinction  conditions 
( aq  =  600  s_1).  The  HCP  temperature  is  varied  in  the  range  1600  <  Thcp  <  1800  K, 
leading  to  super-adiabatic  conditions  with  heat  transfer  from  the  HCP  to  the  reaction  zone. 
The  lower  limit  corresponds  approximately  to  the  extinction  temperature  and  the  upper 
limit  is  just  above  the  adiabatic  flame  temperature,  i.e.  Tq  <  THCp  <  Tad.  The  ignition 
delay  time  was  estimated  to  correspond  to  the  integral  time  scale  of  turbulence  at  a  tem¬ 
perature  of  1200  K  (see  Section  3.6).  For  the  current  Ret,  HCP  supported  re -ignition  is 
therefore  assured. 

3.5.1  Multi-Fluid  Probability  Statistics 

The  multi-fluid  analysis  has  shown  a  clear  redistribution  of  the  fluid  states  with  varying 
Dai.  At  low  Dai,  supported  burning  was  favoured,  while  at  high  Dai  self-sustained  com¬ 
bustion  was  increasingly  realised.  The  reactant  fluid  probability  is  essentially  independent 
of  the  HCP  enthalpy.  The  HCP  impact  becomes  more  evident  for  mixing  fluid  peak  prob¬ 
ability,  see  Fig.  25b.  The  probability  is  reduced  away  from  the  origin  due  to  the  promoted 
reaction  onset  and  significantly  reduced  auto-ignition  delay  time  from  17.6  to  5.1  ps  with 
increasing  Thcp-  This  further  result  in  a  pronounced  augmentation  of  the  mildly  reacting 
fluid  from  13  -  35  %  with  increasing  Thcp  as  depicted  in  Fig.  25c. 

Modest  HCP  blending  fractions  with  insufficient  enthalpy  addition  to  initialise  auto¬ 
ignition  may  cause  elevated  Dai  and  thus  promote  conventional  burning.  Coriton  et  al.  [9] 
observed  the  extension  of  the  intermittent  and  stable  flame  region  towards  higher  product 
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temperatures  for  elevated  a&  of  premixed  CH4  flames.  Although  the  unstable  regime  is  not 
observed  for  the  DME/air  flames  for  the  current  Thcp  range  and  Ret,  a  change  in  burning 
mode  under  super-adiabatic  conditions  is  prominent.  The  augmentation  of  conventional 
burning  is  illustrated  by  means  of  the  strongly  reacting  fluid  probability  in  Fig.  25d,  which 
increases  from  2  %  for  Thcp  =  1600  K  to  approximately  30  %  for  Thcp  =  1800  K.  The 
spatial  extent  of  all  fluid  probabilities  is  essentially  limited  to  one  Lj  for  all  Thcp- 

The  interfaces  encountered  by  traversing  along  the  burner  axis  through  the  quinary 
multi-fluid  field  are  shown  in  Table  7  for  all  cases  investigated  and  visualised  for  Thcp 
=  1600  and  1800  K  in  Fig.  26.  In  the  instance  of  Thcp  =  1800  K,  46  %  of  the  reactant 
fluid  is  observed  adjacent  to  the  mixing  fluid,  which  corresponds  to  20  %  of  all  interfaces. 
The  mildly  and  strongly  reacting  fluid  is  adjoining  the  reactant  fluid  with  a  likelihood  of 
13  and  19  %,  respectively,  which  allocates  23  %  to  the  product  fluid.  At  a  reduced  Thcp 
(=  1600  K),  the  neighbouring  mixing  fluid  probability  increases  to  78  %  with  adjoining 
mildly  and  strongly  reacting  fluids  reduced  to  6  and  1  %,  respectively,  allocating  15  %  to 
the  product  fluid. 

The  mildly  reacting  fluid  is  primarily  (60  %)  in  the  direct  proximity  of  the  product 
fluid  and  independent  of  Thcp-  For  a  high  HCP  enthalpy,  the  residual  adjacent  fluid  is 
composed  of  15  %  reactant  fluid,  16  %  strongly  reacting  fluid  and  8  %  mixing  fluid.  The 
latter  increases  to  24  %  for  Thcp  =  1600  K  at  the  expense  of  reduced  adjacent  reactant 
(11  %)  and  strongly  reacting  (5  %)  fluids.  The  strongly  reacting  fluid  is  found  primarily 
adjacent  to  the  mixing  (24  %)  and  product  (45  %)  fluid  for  the  case  Thcp  =  1600  K,  but  is 
evenly  distributed  (25  ±  3  %)  across  all  fluid  states  for  high  HCP  enthalpies.  Consequently, 
the  interface  statistics  suggest  a  shift  towards  elevated  blending  fraction  for  reduced  HCP 
enthalpy  streams  in  order  to  initialise  the  onset  of  OH  producing  chemical  reactions. 

3.5.2  Conditional  Velocity  Statistics 

The  subsequent  analysis  focuses  on  the  intermediate  cases  with  Thcp  =  1650,  1700  and 
1750  K  due  to  a  modest  disturbance  of  the  stagnation  plane  by  the  LN  HCP  stream  for  the 
limiting  conditions  with  Thcp  =  1600  and  1800  K.  Such  disturbances  can  be  removed  by 
optimising  the  HCP  stream  conditions  for  alternative  temperature  windows.  Unconditional 
velocity  statistics  for  the  discussed  conditions  were  presented  by  Hampp  [22], 

Multi-fluid  conditional  velocity  statistics  are  used  to  delineate  the  influence  of  the  HCP 
enthalpy  on  the  burning  mode  transition.  The  data  are  evaluated  along  the  burner  axis  and 
aligned  at  the  Mie  scattering  iso-contour  (xs  =  0).  The  conditional  mean  axial  fluid  state 
(FS)  velocities  ( Ufs )  are  normalised  by  the  UN  bulk  velocity  ((/{,).  The  axial  (u'0  FS)  and 
radial  (v'0  FS)  velocity  fluctuations  are  normalised  by  the  respective  reactant  fluid  velocity 
fluctuations,  measured  1  mm  away  from  the  nozzle  exit  (u'r  NE,  v'r  NE).  The  normalisation 
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eliminates  modest  discrepancies  between  the  cases  and  provides  a  direct  measure  of  the 
HCP  enthalpy  impact  on  the  turbulent  flow  field. 

The  conditional  reactant  fluid  velocity  prior  the  reaction  onset  (xs  <  0)  is,  naturally 
for  a  flow  with  constant  Dai  and  Ret,  not  affected  by  the  HCP  enthalpy  emerging  the  LN. 
The  impact  of  the  HCP  enthalpy  on  the  thermally  altered  fluid  material  becomes  evident 
in  the  conditional  mixing  fluid  velocity  (Uotmix/Ub)  through  an  increasingly  HCP  driven 
flow  field  for  reduced  Thcp >  see  Fig-  27.  For  Thcp  =  1750  K,  17o,mix/t4  is  negative  in 
the  direct  proximity  of  the  origin,  i.e.  in  line  with  the  reactant  flow,  while  >  0  for  lower 
Thcp-  This  can  be  attributed  to  the  enhanced  blending  capacity  by  the  premixed  reactants 
prior  the  onset  of  OH  producing  chemical  activity  with  reduced  HCP  enthalpy.  The  higher 
HCP  blending  quantities  also  cause  elevated  axial  and  radial  mixing  fluid  fluctuations,  with 
the  latter  being  distinctly  higher  than  the  corresponding  reactant  fluid  velocity  fluctuations 
in  a  comparable  spatial  location. 

In  the  proximity  of  the  origin,  the  conditional  mildly  reacting  fluid  velocity  (Uorrnid 
/ Ub )  is  independent  of  the  HCP  enthalpy  as  illustrated  in  Fig.  28.  The  impact  becomes 
clear  by  the  separation  of  Uo  mUd/Ub  far  away  from  the  origin  (xs  /  Lj  >  0.5).  The 
eased  gradients  for  higher  Thcp  are  caused  by  a  flow  acceleration  towards  the  stagnation 
plane  due  to  a  more  pronounced  dilatation.  The  latter  can  at  least  be  attributed  to  reduced 
HCP  dilution  levels.  The  increasingly  dilatation  driven  flowfield  further  yields  a  distinct 
reduction  in  axial  and  radial  fluctuations.  The  \J (u'u')mud/ \J {u'u’)riNE  show  a  modest 
reduction  compared  to  the  mixing  fluid  while  the  radial  component  exhibits  elevated  fluc¬ 
tuations.  The  conditional  mean  axial  strongly  reacting  fluid  velocity  is,  by  contrast  to  the 
mildly  reacting  fluid  velocity,  insensitive  to  the  HCP  enthalpy.  This  is  characteristic  for  a 
self-sustained  burning  at  constant  Dai  and  Re,  in  the  absence  of  HCP  dilution. 

3.5.3  Conditional  Strain  along  Material  Surfaces 

The  normal  ('«„)  and  tangential  (at  )  strain  conditioned  upon  the  reactant  fluid  iso-contour  is 
depicted  in  Fig.  29  along  with  the  corresponding  total  rate  of  strain  (ad)  and  vorticity  (oj). 
To  recapitulate,  the  LN  HCP  exerted  a  vanishing  effect  on  the  conditional  reactant  fluid 
velocities.  The  an\ i?,  however,  illustrate  a  modest  shift  of  the  PDF  towards  gradually  re¬ 
duced  mean  compressive  values  (—1570  <  an  <  —1250  s_1  for  1650  <  Thcp  <  1750  K). 
The  slight  detachment  of  the  reactant  fluid  iso-contour  from  the  stagnation  plane  with  in¬ 
creasing  HCP  enthalpy  further  causes  a  reduced  mean  tangential  (1000  <  at  <  800  s-1) 
and  total  contracting  (—730  <  <  —500  s-1)  strain  with  reduced  mean  vorticity  levels 

(2810  <  oj  <  2310  s-1).  Moreover,  the  spread  of  all  PDFs  is  reduced  by  ~10  %  with  in¬ 
creasing  Thcp  suggesting  a  first  thermal  alteration  that  is  increasingly  caused  by  adjacent 
exothermic  reactions.  The  velocity  gradients  conditioned  upon  the  mixing  (an\M,  at |  A/, 
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uj | M )  and  mildly  reacting  fluid  surfaces  (an\MR,  cit\MR,  u\MR)  show  similar  trends 
with  details  reported  by  Hampp  [22] . 

For  Thcp  =  1750  K,  the  reaction  onset  of  the  strongly  reacting  fluid  is  anchored  in 
comparatively  low  compressive  strain  regions  with  a  mean  of  an  =  —920  s-1  compared  to 
an  =  -1250  s”1  for  THCp  =  1650  K.  The  PDF(a„| SR)  at  THCp  =  1750  K  further  shows 
a  distinctly  elevated  skewness  towards  decreased  compressive  strain  compared  to  lower 
HCP  enthalpies  as  well  as  to  fluid  states  with  reduced  reactivity.  A  significant  vorticity 
attenuation  of  25  %  is  observed  for  Thcp  =  1750  K  compared  to  1650  K  and  up  to  15  % 
compared  to  the  reactant  fluid  surface. 

3.6  Combustion  Regime  Transitions 

The  estimated  combustion  modes  based  on  Dai  are  shown  in  Fig.  30  using  a  Borghi  dia¬ 
gram  and  are  listed,  along  with  dissipation  rates,  Kolmogorov  length  and  timescales  and 
Kaj  in  Table  2.  The  dissipation  rate  within  the  reactants  was  also  used  to  estimate  a  total 
(mean)  rate  of  strain  ( clt )  via  the  relationship  of  Kostiuk  et  al.  [57]: 

aT  =  ab  +  ai  =  2  ■  (7) 

The  bulk  ab  «  750  s-1  and  turbulent  strain  a/  «  3200  s-1  lead  to  clt  ~  3950  s-1 
(>  aq  for  $  =  1.0)  within  the  reactants  (see  Table  2).  The  rate  of  strain  analysis  along  the 
material  surfaces,  see  Section  3.3.5,  has  illustrated  the  burning  of  conventional  flames  in 
low  strain  regions  with  the  need  of  thermal  HCP  support  to  sustain  combustion  beyond 
the  conventional  (strained  twin  flame)  extinction  limit,  which  exerts  some  influence  for  all 
mixtures. 

The  blending  of  HCP  with  the  reactant  fluid  will  raise  the  temperature  and  the  resulting 
mixture  may  auto-ignite  depending  on  the  delay  time.  A  second  Damkdhler  number  can 
be  used  to  delineate  the  transition  from  conventional  flame  propagation  to  an  auto-ignition 
mode.  Accordingly,  auto-ignition  delay  times  were  computed  for  relevant  initial  tempera¬ 
tures  (1000  -  1700  K)  and  equivalence  ratios  ($  =  0.20,  0.60  and  1.0)  as  shown  in  Fig.  7. 
The  auto-ignition  Damkohler  number  (Da2)  is  here  related  to  the  degree  of  mixing  inside 
a  gas  pocket  as  shown  in  Eq.  (8). 


Da2  =  —  (8) 

Pign 

The  exponential  temperature  dependency  of  Tign  suggests  that  ignition  will  occur  close  to 
the  peak  temperature  within  a  fluid  pocket.  Consequently,  the  expected  range  is  bounded 
by  the  HCP  temperature  of  1700  K  (Da2  =  214;  Tign  =  12.6  pis)  and  the  temperature 
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(«  1196  K;  Tign  =  2.5  ms)  giving  Daa  =  1.  A  corresponding  auto-ignition  Karlovitz 
number  can  readily  be  defined. 


Ka2  =  —  (9) 

Tg 

The  auto-ignition  process  also  depends  on  the  bulk  flow  motion  and  a  third  Damkohler 
number  (Da3),  see  Eq.  (10),  was  defined  based  on  a  convection  residence  time  t/,  (~ 
8.7  ms).  The  latter  was  estimated  by  following  the  mean  trajectory  of  a  fluid  parcel  along 
the  stagnation  point  streamline  from  the  onset  of  reaction/mixing  using  a  reaction  progress 
variable  iso-contour  of  c  =  0.02  [32,  58]  until  it  is  convected  radially  out  of  the  domain  (i.e. 
y=  ±0.015  m). 


Da3  =  —  (10) 

Tign 

The  auto-ignition  Damkohler  numbers  (Da2  and  Da3)  and  the  corresponding  Karlovitz 
number  (Kaa)  were  used  to  derive  the  auto-ignition  regime  diagram  shown  in  Fig.  31.  On 
the  right  hand  side  of  the  Da3  line,  unreacted  mixture  is  likely  to  be  convected  out  of  the 
domain  without  ignition  (assuming  the  absence  of  conventional  flame  propagation).  Left 
of  the  1/Ka2  line,  the  mixture  is  auto-igniting  on  a  timescale  shorter  than  the  Kolmogorov 
timescale.  The  Da3  and  1/Ka2  consequently  bound  the  auto-ignition  manifold  for  present 
study.  The  initial  temperature  (To)  axis  of  this  diagram  can  be  considered  a  third  dimension 
of  a  revised  combustion  regime  diagram  (see  Figs.  30  and  31)  which  intersects  at  a  given 
Re#  and  temperature.  The  conventional  Dai  numbers  can  readily  be  added  to  Fig.  31  at  To 
=  320  K.  Further,  the  impact  of  a  $  variation  on  Da2  is  comparatively  modest  for  DME. 
Under  the  current  conditions  (e.g.  Da3  >  700  and  Ka2  ~  0.05)  any  residual  reactants  will 
auto-ignite  in  the  hot  combustion  products.  It  is  estimated  that  DME  combustion  in  an 
auto-ignition  mode  can  be  sustained  down  to  temperatures  around  1200  K  at  the  current 
Ref.  Further,  no  intermittent  (pulsating)  regime  was  observed.  By  contrast,  pulsating  auto¬ 
ignition  events  have  been  observed  for  CH^air  at  approximately  $  =  0.7  and  Thcp  = 
1700  K,  similar  to  the  observations  by  Coriton  et  al.  [9]. 

The  characteristic  values  for  the  Thcp  variation  are  listed  in  Table  4.  The  conventional 
Dax  is  maintained  constant  at  ~  1.1  with  a  constant  at  Ret  ~  395  and  fine  scale  Karlovitz 
number  (Kai  ~  7.0).  The  auto-ignition  delay  time  ( Tign )  decreases  from  17.6  to  5.1  ps 
within  the  Thcp  range  from  1600  to  1800  K,  which  results  in  a  threefold  reduction  of 
the  auto-ignition  Damkohler  number  (140  <  Da2  <  500).  The  resulting  auto-ignition  Ka2 
increases  gradually  from  0.015  to  0.057  with  decreasing  HCP  enthalpy.  The  results  indicate 
an  impact  on  the  burning  mode  due  to  the  reduced  auto-ignition  delay  time  for  elevated 
HCP  enthalpies.  With  increased  Thcp ,  the  thermochemical  conversion  shifts  away  from  a 
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strong  HCP  dilution  towards  progressively  self-sustained  burning  that  is  possibly  initialised 
by  small  HCP  blending  fractions. 

3.7  Calculation  Methods  for  Distributed  Combustion  Regimes 

The  application  of  calculation  methods  that  do  not  rely  on  presumed  (e.g.  topological  or 
flamelet)  structures  requires  the  introduction  of  multiple  chemical  time-scales.  The  inclu¬ 
sion  of  the  latter  into  turbulent  flowfield  simulations  is  associated  with  exceptionally  high 
computational  costs.  However,  use  of  the  transported  probability  density  function  (pdf)  ap¬ 
proach  for  the  inclusion  of  the  thermochemistry  allows  the  development  of  highly  efficient 
parallel  algorithms  for  large  stochastic  systems.  To  this  effect,  the  scope  of  parallelising 
the  treatment  of  the  chemistry  in  the  context  of  transported  joint-^c/fmethods  has  also  been 
explored  on  different  platforms  in  order  to  permit  more  effective  computations  featuring 
the  direct  implementation  of  extended  mechanisms  for  practical  fuels.  The  latter  are  of  a 
size  approaching  that  of  mechanisms  based  on  the  rapid  decomposition  hypothesis  leading 
to  the  C1-C2  and/or  C1-C4  foundation  chemistry  [26].  The  work  outlined  therefore  forms 
a  step  towards  establishing  if  mechanisms  of  the  latter  type  can  reproduce  measured  trends 
associated  with  combustion  regime  transitions. 

The  Lagrangian  nature  of  stochastic  particle  calculations,  typically  used  in  the  solution 
method  for  the  transported  joint  -pdf,  potentially  permits  an  efficient  parallelisation.  The 
most  computational  demanding  part  of  the  elliptic  calculations  is  the  associated  Monte 
Carlo  part  used  for  the  calculation  of  the  chemical  source  term.  Apart  from  the  molecular 
mixing  model,  the  rest  of  the  method  for  each  stochastic  particle  can  be  calculated  individ¬ 
ually,  which  allows  high  efficiency.  In  addition,  a  high  number  of  particles  of  unmixed  fuel 
or  oxidiser  are  non-reactive  and  thus  have  a  minimal  computational  cost,  which  leads  to 
the  requirement  of  fine  tuning  the  load  balancing  among  the  processes.  In  terms  of  perfor¬ 
mance  cost,  inter-communications  at  the  beginning  of  each  Monte  Carlo  step  have  to  take 
place  in  order  to  distribute  the  particles  among  the  processes.  At  the  end  of  the  Monte  Carlo 
step,  an  operation  also  has  to  take  place  to  calculate  the  global  statistics.  Mixing  models 
form  part  of  the  calculation  method  and  may  be  considered  local  or  non-local.  In  the  case 
of  non-local  models,  the  interacting  particles  are  chosen  randomly  from  the  composition 
space,  while  in  case  of  the  local  models  the  particles  should  mix  with  other  particles  in  their 
immediate  neighbourhood  in  composition  space.  Therefore,  depending  on  the  localness  of 
the  mixing  model  and  the  implementation  of  the  data/particle  parallelisation  (i.e.  how  the 
data/particles  are  stored  among  the  processes),  this  might  lead  into  further  complications 
(e.g.  extra  inter-process  communications). 

Amdahl’s  Law  [64]  can  be  used  to  calculate  the  theoretical  maximum  speed-up  ob- 
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tained  by  using  parallel  processes  as  compared  to  serial  execution.  The  speed-up  can  be 
defined  as  a  relative  change  (or  efficiency),  with  a  potential  maximum  value  of  1  (100  %), 
A  code  can  be  separated  into  two  parts,  the  part  which  can  be  parallelised  and  the  part 
which  cannot  (or  will  not)  be  parallelised.  Let  the  parallelised  part  consume  a  fraction 
”p”  of  the  total  initial  wall  time.  Additionally,  let  ”s”  be  the  reduction  in  latency  of  the 
parallelised  part,  which  in  this  case  is  the  total  number  of  processes.  Then,  the  theoretical 
speed-up  in  term  of  wall  time  is  Si  or  as  a  relative  change  S2- 

Si(s)  =  t - 1 — P  (11) 

1~P~a 

S2(s)=p(l-^J  (12) 

Figure  32  presents  graphically  the  theoretical  maximum  speed-up  as  a  relative  change 
(Eq.  12)  if  the  parallelised  part  of  the  code  accounts  for  50  %,  75  %,  95  %,  and  99  % 
of  the  initial  wall  time.  However,  this  analysis  does  not  take  into  account  various  fac¬ 
tors  which  might  lead  into  increased  latency.  For  example,  using  a  distributed-memory 
communication  implementation,  such  as  MPI,  implies  additional  delays. 

Both  MPI  and  GPGPU  implementations  constitute  Multiple  Data  Stream  Models  [65]. 
Therefore,  how  data  is  stored  in  memory  is  significant  for  their  communications.  Gener¬ 
ally,  there  are  two  distinct  ways  of  storing  the  data:  Structure  of  Arrays  (SoA)  and  Array  of 
Structures  (AoS).  Figure  33  presents  three  actual  code  arrays  OMP,  WTP,  and  PAR  of  size 
250,000  stored  in  a  SoA  on  the  left  and  in  an  AoS  on  the  right.  In  FORTRAN  77,  array  sizes 
are  predefined  at  compilation  time,  therefore  there  is  a  tendency  to  define  significantly  big¬ 
ger  arrays  than  the  necessary.  Thus,  in  blue  colour,  up  to  element  24,000  (arbitrarily  chosen 
as  an  example)  are  the  elements  in  use,  while,  in  orange  colour,  from  element  24,000  until 
the  end  are  noted  as  unused  elements.  The  two  implementations  have  several  trade-offs.  A 
SoA  could  potentially  be  more  effective  since  grouping  together  the  same  database  could 
lead  into  vectorisation.  An  AoS  offers  memory  read/write  coalescing  (all  data  of  an  object 
accessed  at  the  same  time).  However,  if  multiple  processes  trying  to  access  the  same  data 
simultaneously,  then  this  might  lead  into  extra  latency.  Considering  the  data  communi¬ 
cations,  the  AoS  allows  the  bundling  of  all  arrays  into  one  inter-communication  for  each 
process  without  sending  the  unused  elements  of  the  arrays,  the  elements  24,001  to  250,000 
in  Fig.  33.  In  the  current  work,  the  simulation  code  was  rewritten  in  AoS.  The  serial  code 
featuring  SoA  was  roughly  1.5%  slower.  However,  after  its  implementation  on  the  par¬ 
allel  MPI  code,  the  AoS  structure  significantly  reduced  the  MPI  communications  latency, 
improving  remarkably  the  overall  efficiency  as  discussed  below. 
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3.7.1  Efficiency  Gains  from  Parallelisation 

The  current  evaluation  configuration  uses  a  Cabra  burner  geometry  [66,  67]  that  permits 
detailed  studies  of  turbulence-chemistry  interactions.  The  configuration  is  of  particular  rel¬ 
evance  to  auto-ignition  dominated  flows  of  the  kind  that  can  be  expected  to  exert  increasing 
influence  for  distributed  combustion  modes.  Several  studies  have  been  performed  for  this 
geometry  (e.g.  [68,  71])  using  methane  as  fuel.  In  addition  H2  flames  have  been  studied 
and  simulated  using  transported-pc/f  based  methods  [66, 69].  A  wide  range  of  mixing  mod¬ 
els  and  chemical  kinetic  mechanisms  have  also  been  used  with  a  good  agreement  between 
computational  and  experimental  results,  which  suggests  that  the  flames  are  controlled  by 
chemical  kinetics.  There  are,  accordingly,  good  reasons  to  use  this  configuration.  Results 
obtained  with  detailed  methane  chemistry  [71],  featuring  44  species  and  256  reactions,  as 
well  as  a  much  smaller  hydrogen  mechanism  featuring  9  species  and  19  chemical  reactions, 
are  presented  below. 

The  serial  code  version  of  PIPER  ([71])  was  used  as  a  test  bed.  Computer  science 
evolves  rapidly  and  there  have  been  many  changes  regarding  the  Operating  Systems  (OS), 
the  compilers  and  their  libraries.  The  stochastic  nature  of  the  solution  method  makes  it 
potentially  sensitive  to  such  changes.  An  interesting  example  is  the  case  of  extended  pre¬ 
cision,  when  a  single-precision  constant  is  assigned  to  a  double -precision  variable.  It  used 
to  be  evaluated  as  double  precision  according  to  FORTRAN  77  Standards.  However,  ac¬ 
cording  to  Fortran  90/95  Standard  rules,  it  is  now  evaluated  as  single  precision,  which 
essentially  means  that  the  rest  of  the  digits  are  random.  In  order  to  allow  compatibility 
with  FORTRAN  77  Standards,  all  constant  assignments  have  to  be  explicitly  defined  with 
double  precision  accuracy,  or,  in  case  of  Intel  compiler,  the  flag  -fpconstant  can  be  used. 
The  current  testing  was  performed  predominantly  on  the  Imperial  College  HCP  cluster 
operating  under  Oracle  Server  6.7.  However,  other  OS  (OS  X  10.11,  Windows  8.1  and 
CentOS  7)  were  also  evaluated  with  similar  outcomes. 

The  results  from  the  parallelised  code  show,  as  expected,  no  significant  differences 
in  results  as  can  be  seen  in  Fig.  34,  which  compares  the  maximum  OH  concentrations 
of  the  serial  and  the  MPI  versions  of  PIPER  using  the  detailed  mechanisms  for  H2  and 
CH4.  The  differences  between  the  results  are  insignificant  and  can  be  attributed  to  the 
different  Random  Number  Generators  (RNGs).  The  benefits  of  parallelisation  can  be  seen 
in  Fig.  35,  which  presents  the  achieved  speed-up  with  each  bar  representing  the  normalised 
wall  time.  The  run  time  was  hence  reduced  down  to  4.5  minutes  from  27  minutes  in  case 
of  H2  and  down  to  80  minutes  from  nearly  9  hours  in  case  of  CH4.  The  difference  between 
the  theoretical  maximum  speed  and  that  achieved  is  less  than  10%  for  both  cases.  The 
limitations  on  the  number  of  processes  used  is  due  to  the  comparatively  small  physical 
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domains  chosen  (ie.  computational  cells).  Overall,  it  is  clear  that  close  to  the  theoretical 
maximum  speed-up  can  be  achieved.  Accordingly,  it  can  be  expected  that  methods  of  the 
current  type  can  be  combined  directly  with  mechanisms  based  on  the  rapid  decomposition 
hypothesis  leading  to  the  C1-C2  and/or  C1-C4  foundation  chemistry  [26]  to  compute  the 
current  compact  experimental  domains. 

4  Conclusions 

Combustion  regime  transitions  from  close  to  the  corrugated  flamelet  regime  to  distributed 
reaction  zones  were  systematically  investigated  with  examples  given  for  lean  premixed 
DM/air  flames  under  turbulent  conditions  (Ret  ~  380).  An  opposed  jet  back-to-burnt  con¬ 
figuration  with  fractal  grid  generated  turbulence  was  used  with  a  temperature  controlled 
burnt  gas  state.  The  chemical  timescale  was  varied  by  means  of  the  upper  nozzle  equiva¬ 
lence  ratio  (0.2  <  <b[/;Y  <  1.0)  resulting  in  an  estimated  Damkohler  number  (Dai)  range 
of  0.06  <  Da[  <  5.2.  The  stabilising  hot  combustion  products  (HCP)  emerging  from  the 
lower  nozzle  were  created  using  stoichiometric  H2  flames  with  the  temperature  controlled 
by  means  of  CO2  dilution  in  the  range  1600  <  Thcp  <  1800  K. 

The  multi-fluid  concept  of  Spalding  [25]  was  explored  by  means  of  simultaneous  Mie 
scattering,  OH-PLIF  and  PIV.  The  multi-fluid  analysis  has  the  potential  to  extend  methods 
based  on  bimodal  approximations,  such  as  the  BML  [21]  framework,  by  permitting  inter¬ 
mediate  fluid  states,  which  is  of  particular  importance  for  low  Dai  flows.  Five  fluid  states 
were  selected  in  the  present  work  and  thermochemical  interpretations  provided:  reactants, 
combustion  products,  mixing  fluids,  mildly  and  strongly  reacting  (flamelet-like)  fluids. 

The  combustion  regime  transition  from  self-sustained  flame  propagation  towards  a 
distributed/supported  burning  mode  was  quantified  using  the  multi-fluid  analysis  with  a 
gradual,  but  distinct,  redistribution  between  the  defined  fluid  states  observed  with  decreas¬ 
ing  Dai.  The  interface  and  adjacent  fluid  state  statistics  provided  further  insight  of  the 
observed  fluid  state  reallocation.  The  fluid  state  that  represents  self-sustained  burning  was 
observed  with  a  40  %  probability  directly  adjacent  to  the  reactant  fluid  for  $  =  1.0  (Dai  = 
5.1),  which  was  significantly  reduced  with  Dai  in  favour  of  interlayers  of  reduced  reactiv¬ 
ity.  The  less  reactive  fluid  states  were  increasingly  separated  from  the  reactants  by  a  layer 
acting  as  an  enthalpy  source.  This  resulted  in  a  mildly  reacting  fluid  peak  probability  that 
was  shifted  towards  the  stagnation  plane. 

The  current  work  further  showed  a  distinct  fluid  acceleration  due  to  increased  heat 
release  at  high  Dai,  while  the  flow  field  was  dominated  by  the  hot  combustion  product 
blending  for  less  reactive  mixtures.  The  multi-fluid  delineation  of  the  combustion  process 
hence  allowed  the  separation  of  the  contributions  made  by  the  various  fluid  states.  The 
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proximity  of  the  stagnation  plane  was  characterised  by  a  high  compressive  axial  strain 
that  increasingly  prevented  self-sustained  burning  for  reduced  Damkohler  numbers.  The 
increasing  influence  of  hot  combustion  product  support  resulted  in  highly  strained  and  Dai 
independent  chemically  active  zones  with  no  favoured  dilatation  direction  -  consistent  with 
an  auto-ignition  related  burning  mode.  By  contrast,  a  preferential  alignment  of  the  flame 
front  normal  with  the  principal  extensive  strain  and  significantly  reduced  vorticity  levels 
were  evident  for  self-sustained  flames  with  high  Dai  ■  The  associated  dilatation  caused  a 
reduction  of  the  gradient  scalar  flux.  However,  the  transition  to  counter-gradient  scalar 
transport  was  suppressed  under  the  current  conditions  (Re4  ~  380). 

The  multi-fluid  analysis  was  further  used  to  delineate  the  impact  of  the  HCP  enthalpy 
on  the  burning  mode  transition.  The  DME  case  marking  the  estimated  transition  from  the 
thin  reaction  zone  to  distributed  burning  (i.e.  <f>  =  0.6)  was  stabilised  against  HCP  streams 
of  varying  temperature  (1600  <  Thcp  <  1800  K).  The  multi-fluid  probability  and  inter¬ 
face  statistics  illustrated  the  advanced  initialisation  of  chemically  active  fluid  states  with 
reduced  blending  fractions  for  higher  Thcp-  A  substantial  augmentation  of  the  strongly 
reacting  fluid  was  observed  for  Thcp  in  excess  of  the  adiabatic  flame  temperature.  By 
contrast,  the  strongly  reacting  fluid  nearly  vanished  for  Thcp  of  the  order  of  the  extinc¬ 
tion  temperature.  The  absence/reduction  of  HCP  diluents  in  the  chemically  active  fluids 
caused  a  transition  away  from  distributed  towards  self-sustained  burning  with  increasing 
HCP  enthalpy.  The  turbulence  -  chemistry  interaction  was  evaluated  in  the  range  1650  < 
Thcp  <  1750  K.  Despite  the  narrow  variation,  the  exponential  temperature  dependency  of 
the  auto-ignition  process  resulted  in  clear  trends  and  highlighted  an  incrementally  stronger 
and  directed  dilatation  with  increasing  HCP  enthalpy  that  was  attributed  to  reduced  blend¬ 
ing  fractions.  By  contrast,  the  velocity  statistics  of  the  strongly  reacting  fluid,  representing 
a  fluid  state  with  little  or  no  HCP  support,  was  not  affected  by  the  HCP  enthalpy.  The  strain 
distribution  along  the  multi-fluid  material  surfaces  indicated  a  gradual  detachment  of  the 
reaction  onset  from  the  stagnation  plane  with  increasing  HCP  enthalpy  that  was  attributed 
to  an  eased  reaction  onset  at  lower  HCP  blending  fractions. 

The  results  of  all  cases  were  further  evaluated  in  the  context  of  conventional  and  auto¬ 
ignition  based  Damkohler  and  Karlovitz  numbers  and  it  was  shown  that  the  failure  to  estab¬ 
lish  conventional  turbulent  flames  leads  to  a  combustion  mode  that  is  dominated  by  the  hot 
combustion  product  temperature.  With  increasing  HCP  enthalpy,  self-sustained  burning 
was  promoted  due  to  an  eased  reaction  onset.  The  classification  of  the  data  in  a  conven¬ 
tional  Borghi  diagram  showed  the  wide  distribution  with  varying  equivalence  ratio,  while 
the  auto-ignition  dimensions  advanced  the  fundamental  understanding  of  the  observed  dif¬ 
ferences  of  the  turbulence-chemistry  interaction  between  the  investigated  fuels. 

The  current  work  has  illustrated  the  need  for  a  more  comprehensive  analysis  technique 
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for  combustion  processes  with  low  Damkohler  and  high  turbulent  Reynolds  number.  The 
current  multi-fluid  delineation  removes  the  assumption  of  thin  reaction  zones  and  has  pro¬ 
vided  novel  insight  into  combustion  regime  transitions.  The  analysis  further  highlighted 
the  distinct  need  to  accurately  model  the  combustion  chemistry  -  in  particular  for  low  Dai 
flows.  The  obtained  data  are  expected  to  be  particularly  useful  for  the  evaluation  of  models 
that  do  not  assume  a  flame  front  structure  and  may  thus  be  used  to  characterise  burning 
mode  transitions.  It  has  been  shown  that  such  models  can  be  efficiently  parallelised  with  a 
performance  close  to  the  theoretical  maximum. 
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Tab.  1:  Summary  of  utilised  OH  intensities  (7 ..)  for  varying  fuel  and  equivalence  ratio  to 
define  the  thresholds  A oH,prod,<s>  and  A oH,str,<s>- 
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Fluid  state  normalised  signal  intensities. 

$ 

0.2 

0.4 

0.6 

0.8 

1.0 

P 

1.0 

Ilimit 

<0.5 

Ihcp 

<  2.0 

Ibcp,<s> 

<  2.0 

3.0 

3.6 

Iprod^ 

1.0 

-2.0 

1.0 -3.0 

1.0 -3.6 

>2.0 

2.0 -3.4 

3.0 -5.6  3.6 -6.2 

Istr,& 

N.A. 

>  3.4 

>  5.6 

>  6.2 

A.OH,prod,& 

2.0 

2.0 

2.0 

3.0 

3.6 

A.QH,str,& 

N.A. 

N.A. 

3.4 

5.6 

6.2 

Tab.  2:  Summary  of  turbulent  and  chemical  conditions  to  derive  the  turbulent  Reynolds, 
Damkohler  and  Karlovitz  numbers  for  DME/air  mixtures  at  varying  <f>  and  low 
rate  of  strain  (a  =  75  s’  1 ).  The  turbulence  conditions  are  evaluated  within  the 
reactants.  The  auto-ignition  delay  time  ( rign ),  Da2,  Ka2  and  Da3  are  evaluated  at 
1700  K. 
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2.46 

ms-1 

0.25 

0.24 

0.23 

0.22 

0.22 

Lr/ 

jjm 

71 

71 

73 

75 

74 

TV 

ms 

0.29 

0.29 

0.31 

0.33 

0.33 

Tb 

ms 

8.68 

aq 

S"1 

N.A. 

N.A. 

600 

2000 

3100 

clt 

s-1 

4200 

4160 

3940 

3750 

3755 

£r 

m2  s-3 

210 

200 

173 

150 

151 

ry  X  106 

m2  s_1 

17.5 

17.2 

17.0 

16.7 

16.5 

Re, 

- 

373 

379 

383 

378 

415 

Dai 

- 

0.08 

0.14 

1.18 

3.86 

5.13 

Da2 

- 

214  ±  13 

Da3 

- 

701  ±  42 

Ka! 

- 

112 

62.8 

7.06 

2.09 

1.45 

Ka2 

- 

0.046  ±  0.003 
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Tab.  3:  Summary  of  utilised  OH  intensities  (I...)  for  DME/air  at  $  =  0.6  at  varying  Thcp 
to  define  the  thresholds  A oH,prod,T  and  A oH,str,T- 


Fluid  State 

Fluid  state  intensities  . 

Thcp  [K] 

1600 

1650 

1700 

1750 

1800 

Imild,T' 

2.0 -3.9 

2.0 -3.6 

2.0 -3.3 

2.0 -3.0 

2.0 -2.8 

Istr,T 

>  3.9 

>  3.6 

>  3.3 

>  3.0 

>  2.8 

^OH,Prod,T 

2.0 

2.0 

2.0 

2.0 

2.0 

J^OH,str,T 

3.9 

3.6 

3.3 

3.0 

2.8 

Tab.  4:  Summary  of  turbulent  and  chemical  conditions  to  derive  the  turbulent  Reynolds, 
Damkohler  and  Karlovitz  numbers  for  DME/air  at  4>  =  0.6  at  low  strain  rate  (a  = 
75  s-1).  The  turbulence  conditions  are  evaluated  within  the  reactants.  The  auto¬ 
ignition  delay  time  ( Tign ),  Da2,  Ka2  and  Da3  are  evaluated  at  the  respective  Thcp- 
The  data  are  supplementary  to  the  data  in  Table  2. 


Thcp 

K 

1800 

1750 

1700 

1650 

1600 

1~ign 

ps 

5.07 

6.66 

8.96 

12.4 

17.6 

U rms 

ms-1 

1.63 

1.67 

1.65 

1.60 

1.64 

£r 

m2  s-3 

158 

178 

172 

183 

176 

Ret 

- 

393 

403 

397 

385 

395 

Dai 

- 

1.15 

1.12 

1.13 

1.17 

1.14 

Da2 

- 

495 

369 

277 

206 

142 

Da3 

- 

1712 

1303 

969 

700 

493 

Kai 

- 

6.67 

7.09 

6.96 

7.19 

7.05 

Ka2 

- 

0.015 

0.021 

0.029 

0.041 

0.057 

Tab.  5:  Physical  and  resolved  length  scales  in  pm.  A piy  and  A m f  is  the  PIV  and  multi¬ 
fluid  resolution,  respectively.  The  current  minimum  laminar  flame  thickness  at  4> 
=  1.0  is  mini  Sf).  A  d  is  the  mean  scalar  dissipation  layer  thickness  [59]  and  A  b  the 
Batchelor  scale  [60]. 


Scales 

Reactants 

HCP 

A  piv 

300 

Am  f 

250 

min  (5f) 

222 

N.A. 

Ad 

621  ±  18 

N.A. 

A  B 

86  ±  15 

~  370 
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Tab.  6:  Multi-fluid  interface  probabilities  for  DME/air  mixtures  in  percent  by  traversing 
along  the  stagnation  point  streamline  in  3000  instantaneous  images  from  reactants 
(R)  to  products  (P).  M  -  mixing  fluid,  MR  -  mildly  reacting  fluid,  SR  -  strongly 
reacting  fluid.  Interfaces  marked  with  (-)  are  not  defined. 


Interface 

0.2 

0.4 

$ 

0.8 

1.0 

0.6 

Samples  xlO-3 

7.5 

7.8 

oo 

oo 

9.4 

9.8 

R 

M 

43 

38 

36 

10 

5 

R 

MR 

<1 

<1 

<1 

<1 

2 

R 

SR 

- 

- 

1 

10 

14 

R 

P 

2 

2 

1 

12 

12 

M 

MR 

1 

6 

8 

<1 

<1 

M 

SR 

- 

- 

10 

2 

1 

M 

P 

39 

30 

15 

8 

5 

MR 

SR 

- 

- 

7 

15 

18 

MR 

P 

15 

24 

21 

25 

27 

SR 

P 

- 

- 

3 

17 

15 

Tab.  7 :  Multi-fluid  interface  probabilities  for  varying  HCP  enthalpy  in  percent  by  travers¬ 
ing  along  the  stagnation  point  streamline  in  3000  instantaneous  images  from  re¬ 
actants  (R)  to  products  (P).  M  -  mixing  fluid,  MR  -  mildly  reacting  fluid,  SR  - 
strongly  reacting  fluid.  Interfaces  marked  with  (-)  are  not  defined. 


Interface 

1800 

1750 

Thcp  [K] 

1700 

1650 

1600 

Samples 

8953 

10018 

9489 

8349 

8390 

R 

M 

20 

27 

27 

31 

36 

R 

MR 

6 

3 

4 

4 

3 

R 

SR 

8 

3 

3 

2 

<1 

R 

P 

10 

6 

6 

6 

7 

M 

MR 

3 

4 

7 

9 

7 

M 

SR 

7 

8 

7 

4 

1 

M 

P 

10 

15 

13 

18 

28 

MR 

SR 

6 

6 

6 

4 

1 

MR 

P 

22 

22 

24 

22 

16 

SR 

P 

8 

4 

2 

1 

<1 

DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


8  Figures 


46 


8  Figures 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


8  Figures 


47 


Fig.  1 :  Experimental  configuration. 


Unburnt  Reactants 

t* 


Hot  Combustion  Products 
(HCP) 


Fig.  2:  Coordinate  system  convention  and  location  of  reference  windows  (dashed  rectan¬ 
gles). 
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Fig.  3:  Identification  of  the  multiple  fluid  states:  (a)  Mie  scattering;  (b)  Planar  OH  LIF 
image.  PIV  vectors  overlaid.  Pink  line  is  the  DS  iso-contour  enclosing  the  reactant 
fluid,  yellow  line  encloses  the  entire  OH  field,  green  line  encloses  the  mildly  react¬ 
ing  fluid,  and  red  line  encloses  the  strongly  reacting  (flamelet)  fluid.  The  mixing 
fluid  is  bounded  between  the  pink  and  yellow  line. 
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Fig.  4:  Schematic  visualisation  of  the  multi-fluid  post-processing  methodology.  PSD  - 
particle  seeding  density;  I  -  experimental  OH  -  PLIF  signal  intensity. 
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Fig.  5:  Perfectly  stirred  reaction  calculations  to  evaluate  the  OH  concentration  following 
mixing  of  cold  air  with  HCP  at  varying  Thc p  for  the  integral  (Yj)  and  Kolmogorov 
(t„)  timescale  to  establish  the  maximum  OH  signal  level  in  the  absence  of  fuel. 
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Fig.  6:  Theoretical  mixture  seeding  density  due  to  UN  and  LN  stteam  mixing  as  a  function 
of  blending  quantity  (i.e.  mixing  temperature)  and  seeding  density  ratio  for  Thcp 


=  1700  K. 


To  [K] 


1/T0  •  103  [K-1] 


Fig.  7:  Perfectly  stirred  reaction  calculations  to  evaluate  auto-ignition  delay  times  for 
DME-air  with  the  bar  representing  the  variation  with  $  =  0.2  -  1 .0,  P  =  1  atm 
and  validation  calculations  at  $  =  1.0,  P  =  1.0  -  1.8  bar  with  experimental  auto¬ 
ignition  delay  times  (□  -  [61]  $  =  1.0,  P  =  1.8  atm;  O  ~  [62]  $  =  1.0,  P  =  1.0  bar; 
x  -  [63]  $  =  1 .0,  P  =  1.2  atm)  for  1  %  DME  in  CU/Ar  mixtures. 
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Fig.  8:  Laminar  flame  calculations  in  a  twin  (left  column)  and  back-to-burnt  (right  col¬ 
umn)  configuration  to  determine  the  extinction  strain  (vertical  dashed  arrows)  and 
to  define  the  strongly  reacting  fluid  threshold  A oH,str,<s>  (horizontal  dashed  ar¬ 
rows):  1st  row:  Maximum  temperature;  2nd:  Peak  OH  concentration. 
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Strain  Rate  [s  ^  Strain  Rate  [s  ^ 

Fig.  9:  Laminar  flame  calculations  to  visualise  the  change  in  burning  mode  with  increasing 
strain  in  the  back-to-burnt  configuration  for  premixed  DME/air  flames  at  varying 
<f>.  Top:  Flame  thickness  Sf  based  on  the  5  -  95  %  of  the  fuel  consumption; 
Bottom:  Laminar  burning  velocity. 
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Fig.  10:  Laminar  flame  calculations  in  a  BTB  geometry  for  varying  HCP  conditions,  i.e 
LHS  BC:  HCP  with  Thcp  and  respective  species  concentrations;  Constant  RHS: 
premixed  DME/air  at  $  =  0.6  (Dai  =  1.0).  (a)  Laminar  flame  thickness  based 
on  the  95  -  5  %  fuel  consumption;  (b)  Laminar  burning  velocity;  (c)  Normalised 
OH  peak  concentration;  (d)  Damkohler  number. 
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Fig.  1 1 :  OH  lifetime  (t)  in  the  post  combustion  environment  for  $  =  0.60  and  1 .00  under 
extinction  conditions.  The  OH  radical  drops  within  ~3.0  r^  to  the  OH  level  of 
the  HOP  fluid  f[OH]ffCp  =  2.0 -OH*)  for  the  case  $  =  0.60  and  ~3.6  OH*  for 
$  =  1.00  respectively.  The  product  fluid  thresholds  (A oh, prod.,®)  are  marked  and 
A OH,str,®  can  be  inferred  at  t  /  tj  =  0. 


Fig.  12:  Sensitivity  analysis  on  thresholds  definition:  (a)  Product  fluid  threshold 
(A oH,prod,^>)  and  (b)  strongly  reacting  (flamelet)  fluid  threshold  (A oH,str,®)  for 
$  =  0.6.  Arrows  indicate  increasing  threshold  value. 
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Fig.  13:  Example  of  a  quinary  multi-fluid  field  (DME/air- $  =  0.8):  Light  blue  -  reactant 
fluid,  dark  blue  -  mixing  fluid,  orange  -  mildly  reacting  fluid,  red  -  strongly 
reacting  fluid  and  green  -  product  fluid.  The  pink  iso-contour  encloses  the  re¬ 
actant  fluid,  the  white  iso-contour  encloses  all  OH  signal,  the  black  iso-contour 
encloses  the  mildly  reacting  fluid,  the  red  iso-contour  encloses  the  strongly  re¬ 
acting  (flamelet)  fluid  and  the  mixing  fluid  is  bordered  by  the  pink  and  white 
iso-contour. 
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Fig.  14:  Multi-fluid  probability  statistics  for  DME  cases  at  $  =  0.2  -  1.0  and  aligned  at 
the  first  thermal  alteration  iso-contour:  (a)  Reactant  fluid;  (b)  Mixing  fluid;  (c) 
Mildly  reacting  fluid;  (d)  Strongly  reacting  (flamelet-like)  fluid  probability.  The 
markers  on  the  DME  -  <t>  =  0.20  line  are  drawn  for  identification  purposes  and  do 
not  represent  the  actual  resolution. 
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Fig.  15:  Multi-fluid  interface  diagram  for  DME/air  mixtures  at  $  =  0.2  (left),  <f>  =  0.6 
(middle)  and  $  =  1.0  (right)  and  elaborated  for  all  cases  in  Table  6.  The  distinct 
shift  towards  higher  reactive  and  more  exothermic  intermediate  fluid  states  is 
illustrated  by  means  of  the  probability  weighted  arrows,  indicating  adjacent  fluid 
states. 
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Fig.  16:  Unconditional  velocity  and  its  fluctuation  along  the  stagnation  point  streamline: 

Mean  axial  velocity  (top  row),  axial  velocity  fluctuation  (middle  row)  and  ra¬ 
dial  velocity  fluctuation  (bottom  row).  Only  every  second  data  point  is  drawn  to 
enhance  the  readability. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


8  Figures 


60 


-0.25  0  0.25  0.5 

X  /  D  [-] 

Fig.  17:  Conditional  mean  axial  reactant  fluid  velocity  and  fluctuations  along  the  stag¬ 
nation  point  streamline:  Top  -  Ur/Ub\  Middle  -  yj ( u'u')r/Ub ;  Bottom  - 
\/( v'v')r/Ub . 
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Fig.  18:  Conditional  mean  axial  mixing  fluid  velocity  and  fluctuations  along  the  stagna¬ 
tion  point  streamline:  Top  -  UmiX/Ub',  Middle  -  \J ( u'u')miX/Ub ;  Bottom  - 
\J  ( v r  V  !  )  m  7  x  /  U  b  - 
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Fig.  19:  Conditional  mean  axial  mildly  reacting  fluid  velocity  and  fluctuations  along  the 
stagnation  point  streamline:  Top  -  Urnlid/Ui,;  Middle  -  yj ( u'u')mud/Ub\  Bottom 

“  \/ iy^^mild/Ub. 
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Fig.  21:  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  reactant  fluid  surface: 

Top  left:  Normal  strain;  Right:  Tangential  strain;  Bottom  left:  Total  strain;  Right 
Vorticity.  Legend  refers  to  <h. 


-6  -4  -2  0  2  4  -2  0  2  4  6 


-4  -2  0  2  4  -2  0  2  4  6  8 

ad  x  10‘3  [s'1]  u  x  10'3  [s'1] 


Fig.  22:  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  mildly  reacting  fluid 
surface:  Top  left:  Normal  strain;  Right:  Tangential  strain;  Bottom  left:  Total 
strain;  Right  Vorticity.  Legend  identical  to  Fig.  21. 
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Fig.  23:  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  strongly  reacting  fluid 
surface:  Top  left:  Normal  strain;  Right:  Tangential  strain;  Bottom  left:  Total 
strain;  Right  Vorticity.  Legend  identical  to  Fig.  21. 
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Fig.  24:  Limitation  of  a  bimodal  description  evaluated  by  the  difference  (A)  between  a 
particle  density  and  OH  signal  based  segregation  technique  along  the  reaction 
progress  c.  Discrepancy  of  the  mean  axial  reactant  fluid  (top  left),  product  fluid 
(top  right)  and  slip  (bottom  left)  velocity  as  well  as  scalar  flux  (bottom  right). 
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Fig.  25:  Multi-fluid  probability  for  DME  cases  at  $  =  0.6  with  Thcp  =  1600  -  1800  K: 
(a)  Reactant  fluid;  (b)  Mixing  fluid;  (c)  Mildly  reacting  fluid;  (d)  Strongly  re¬ 
acting  (flamelet-like)  fluid  probability.  The  markers  are  drawn  for  identification 
purposes  and  do  not  represent  the  actual  resolution. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


8  Figures 


67 


Fig.  26:  Multi-fluid  interface  diagram  for  DME  /  air  at  $  =  0.6  and  a  supporting  HCP 
temperature  of:  (a)  Thcp  =  1600  and  (b)  1800  K.  A  shift  towards  higher  HCP 
blending  fraction  for  reduced  HCP  enthalpy  is  required  to  initialise  OH  producing 
chemical  activity  with  absolute  interface  probabilities  listed  in  Table  7. 
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Fig.  27:  Conditional  mean  axial  mixing  fluid  velocity  and  fluctuations  along  the  stag¬ 
nation  point  streamline  and  aligned  at  the  Mie  scattering  iso-contour  (xs 

=  0):  Top  -  Uo,mix/Ub\  Middle  -  sj (u'u')0^mix/ ^ (uV)r, ne\  Bottom  - 

\J  (f,'l’,)o,mix/v/  [v'v')r,NE- 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


8  Figures 


69 


0.5 


0.4 

0.3 

0.2 

0.1 


0 

1 


K) 

I? 

f-T 

>  0.8 


0.4 


^  1.8 

c 

>  1.6 


s> 


1 


0  0.5  1  1.5 


xs/Li[-j 


Fig.  28:  Conditional  mean  axial  mildly  reacting  fluid  velocity  and  fluctuations  along  the 
stagnation  point  streamline  and  aligned  at  the  Mie  scattering  iso-contour  (xs 

=  0):  Top  -  U0,miid/Ub\  Middle  -  sj {u'u'^miid/ \J {u'u')r,NE\  Bottom  - 

\J  Wv')o,miid/  \/(v'v')r}NE- 
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Fig.  29:  PDF  of  the  rate  of  strain  and  vorticity  evaluated  along  the  reactant  fluid  iso¬ 
contour:  Top  left:  Normal  strain;  Right:  Tangential  strain;  Bottom  left:  Total 
strain;  Right  Vorticity. 


Fig.  30:  Combustion  regime  transitions  of  DME/air  mixtures  at  varying  $  visualised  in  a 
Borghi  diagram.  Symbols  as  in  Figs.  16  -  20. 
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Fig.  31:  Auto-ignition  manifold  for  DME/air  mixtures  with  Damkohler  numbers  based 
on  the  turbulent  integral  timescale  (Daa),  the  bulk  flow  timescale  (Da3)  and  the 
inverse  of  the  auto-ignition  Karlovitz  number  (1/Ka2)  as  a  function  of  initial  tem¬ 
perature  (To).  The  conventional  Dai  for  the  different  equivalence  ratios  (symbols 
as  in  Figs.  16  -  20)  are  also  shown  at  T0  =  320  K.  The  bars  on  the  Da2  line 
indicate  the  mixture  impact  ($  =  0.2  -  1.0). 
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Fig.  32:  The  theoretical  maximum  speed-up  according  to  Amdahl’s  Law  [64]. 
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Fig.  33:  The  difference  in  memory  alignment  of  a  Structure  of  Arrays  and  Arrays  of  Struc¬ 
tures. 
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Fig.  34:  Comparing  the  maximum  OH  concentrations  obtained  with  the  serial  version  of 
PIPER  and  that  parallelised  using  MPI.  The  left  figure  was  obtained  with  a  de¬ 
tailed  mechanism  for  a  H2  flame  and  the  right  figure  with  a  detailed  mechanism 
for  the  corresponding  CH  |  flame. 
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Fig.  35:  The  achieved  speed-up  of  PIPER  using  MPI.  The  left  figure  feature  a  detailed 
mechanism  for  H2  and  the  right  figure  the  detailed  mechanism  for  the  corre¬ 
sponding  CH4  flame.  The  tail-off  in  gain  is  due  to  the  small  domain  size  (number 
of  computational  cells  used). 
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